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The  research  that  led  to  the  developments  reported  herein  was  conducted 
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was  also  sponsored  by  the  U.  S.  Department  of  Agriculture  Forest  Service. 
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Messrs.  E.  C.  Odom  and  R.  W.  Patrick,  under  the  general  supervision  of 
Mr.  J.  P.  Sale,  Chief  of  the  Geotechnical  Laboratory,  WES.  This  report  is 
essentially  Dr.  Barber’s  dissertation,  which  was  submitted  to  Texas  A&M 
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Philosophy  degree.  The  authors  were  provided  editorial  assistance  by 
Mrs.  L.  M.  Beall. 

Directors  of  WES  during  the  investigation  were  COL  G.  H.  Hilt,  CE,  and 
COL  J.  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 
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INTRODUCTION 


Background 

Pavement  design  has  traditionally  been  set  apart  from  maintenance 
and  repair  considerations.  Early  design  procedures,  by  omission,  pro- 
vided for  the  design  of  pavements  that  were  to  perform  some  service  up 
to  a point  where  sudden  failure  occurred  as  a result  of  some  predeter- 
mined quantity  of  traffic.  This  "failure"  was  some  finite  definition 
of  the  pavement  condition.  However,  it  has  always  been  intuitively  ob- 
vious that  the  pavement  deterioration,  or  damage,  began  to  occur  upon 
initiation  of  traffic  and  gradually  accrued  to  some  point  where  condi- 
tions were  unsatisfactory.  It  was  also  obvious  to  the  designer  that 
this  unsatisfactory  condition  varied  from  one  facility  to  another  de- 
pending upon  the  needs  and  desires  of  the  user. 

More  recently,  and  especially  in  the  past  decade,  designers  have 
sought  methodology  to  quantify  deterioration  of  pavements  in  various 
modes  and  to  properly  define  failure  of  a pavement.  These  achievements 
have  been  considered  essential  in  order  to  not  only  design  a pavement 
but  obtain  the  highest  possible  benefit  from  a pavement  throughout  its 
entire  life. 

This  concept  has  been  termed  "life-cycle  management"  of  pavements. 
Life-cycle  management  of  pavements  is  a new  concept,  but  it  is  clearly 
based  upon  the  classical  definition  of  engineering  itself.  Life-cycle 
management  can  be  described  as  the  management  of  a pavement  from  its 


The  citations  on  the  following  pages  follow  the  style  of  the  Journal 
of  the  Geotechnical  Engineering  Division,  Proceedings  of  American  Society 
of  Civil  Engineers. 
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inception  until  the  end  of  its  life.  The  term  then  should  include  plan- 
ning, design,  maintenance,  repair,  and  some  control  of  usage. 

In  life-cycle  management,  essentially  every  tradeoff  must  he  opti- 
mized. Particularly,  a design  must  be  at  some  minimum  cost  with 
respect  to  deterioration.  Additionally,  the  design  must  be  aimed  di- 
rectly at  the  level  of  reliability  desired  as  a procedure  for  minimizing 

V 

design  redundancy  or  design  insufficiency. 

The  U.  S.  Army  Corps  of  Engineers  (CE)  has  been  typical  in  their  his- 
tory of  developing  pavement  design  procedures.  During  the  early  19^0 's  the 
CE  adopted  the  California  Bearing  Ratio  (CBR)  tests  for  defining  material 
strength  and  developed  the  CE  design  method  for  flexible  pavements  (1,2).^ 
This  semiemphirical  method  was  selected  in  a time  of  military  need  partly 
for  its  ease  of  application.  Similarly,  a theoretical  analysis  character- 
ized by  an  elastic  plate  on  a liquid  subgrade  was  selected  for  development 
of  the  CE  design  method  for  rigid  pavements.  The  foundation  strength  for 
rigid  pavement  design  was  characterized  by  the  modulus  of  subgrade  reac- 
tion (16) . The  CE  design  methods  for  rigid,  flexible,  and  other  pavement 
types  have  undergone  several  modifications  through  the  years  and  are  still 
in  use  today.  These  methods  are  deterministic  in  that  they  provide  for 
design  of  new  pavements  with  respect  to  a specific  failure  cirterion,  but 
not  for  analysis  of  gradual  deterioration  in  the  respective  modes,  or 
"system  drift."  This  implies  sudden  failure  of  pavements  upon  application 
of  some  computed  quantity  of  traffic  applications. 

Major  research  programs  have  been  undertaken  in  the  ast  decade  to 

^Numerals  in  paraentheses  refer  to  corresponding  items  in  Appendix 
I. — References. 
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improve  CE  pavement  design  and  management  capabilities.  In  the  early 
1970' s the  WES  initiated  research  programs  to  improve  design  procedures 
for  the  CE  and  other  agencies.  Theoretical  design  procedures  using 
basic  material  parameters  have  been  developed  as  a result  of  research 
sponsored  by  the  CE  and  the  Federal  Aviation  Administration  (FAA). 

These  procedures  apply  chiefly  to  rigid  (25)  and  flexible  (h)  airport 
and  road  pavements  as  well  as  some  variations  of  these  pavement  types. 

However,  these  improved  procedures  are  still  deterministic  in  that 
they  utilize  a specific  failure  criterion,  assume  one-time  failure  at 
some  point,  and  do  not  address  rate  of  deterioration  or  pavement 
reliability. 

A program  of  study  to  fulfill  the  needs  of  the  CE  and  FS  was 
approved  and  funded  in  197*+  to  develop  life-cycle  procedures  for 
pavements  based  upon  pavement  deterioration  and  statistical 
reliability  (3,10).  This  research  is  aimed  at  quantification  of  de- 
terioration of  pavements  and  assessing  reliability.  This  program  of 
study,  currently  in  progress,  is  partially  sponsored  by  the  U.  S.  De- 
partment of  Agriculture  Forest  Service  (FS).  The  FS  is  participating 
in  the  study  as  a result  of  their  need  to  assess  damage  to  roads  caused 

j 

by  logging  operations  (3).  Data  available  or  being  collected  are  ex- 
pected  to  provide  for  analysis  of  damage  caused  by  various  types  of 
vehicles.  This  capability  will  provide  the  basis  for  development  of  a 
differential  cost  analysis  procedure  to  aid  in  assessment  of  maintenance 
costs  to  private  sector  timber  industry. 

The  Department  of  Defense,  through  the  CE,  is  participating  in  the 
effort  as  a result  of  determination  of  the  need  not  only  to  predict 

3 
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deterioration  of  roads  in  a military  scenario  but  to  assess  pavement  re- 
liability in  military  tactical  and  logistical  operations.  This  capa- 
bility, along  with  current  reliability  concepts  that  are  being  applied 
to  military  vehicle  operation,  will  provide  for  reliability  of  the  over- 
all systems  that  include  ground  and  air  vehicles  as  well  as  the  mediums 
(pavements)  upon  which  they  operate. 

Objectives  of  This  Study 

The  objective  of  this  study  is  to  investigate  the  hypothesis  that 
effective  pavement  life-cycle  management  can  be  achieved  through  utiliza- 
tion of  deterioration  and  reliability  concepts.  In  order  to  accomplish 
this  investigation,  several  intermediate  objectives  are  set  forth  as 
follows : 

1.  Utilize  the  surface  rutting  mode  of  deterioration  to  develop  a 
pilot  deterioration  prediction  procedure. 

2.  Further  develop  the  rutting  prediction  procedure  into  a relia- 
bility assessment  system. 

3.  Combine  the  deterioration  and  reliability  models  into  a deterio- 
ration and  reliability  analysis  procedure  for  use  in  life-cycle  manage- 
ment. This  procedure  will  be  a pilot  procedure  that  incorporates  the 
rutting  mode  of  deterioration. 

k.  Provide  a basis  for  expanding  these  developments  to  include 
other  modes  of  deterioration. 

Accomplishment  of  these  objectives  will  establish  a basis  for  the 
development  of  an  effective  life-cycle  management  procedure. 
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Scope  of  Work 

The  initial  efforts  consisted  of  a search  for  existing  data  to  de- 
termine whether  enough  rutting  data  were  available  to  provide  a basis 
for  development  of  the  deterioration  and  reliability  models.  The  data 
were  analyzed  and,  being  found  tentatively  satisfactory,  were  utilized 
for  this  purpose. 

Literature  was  reviewed,  and  studies  were  made  to  determine  the 
most  suitable  method  of  analysis  of  the  available  data.  After  selecting 
a method  for  data  analysis,  a major  portion  of  the  research  effort  con- 
sisted of  analysis  of  data  and  comparison  of  existing  data  with  that 
being  accumulated  in  ongoing  field  evaluations. 

As  deterioration  and  reliability  models  came  forth  from  tha  anal- 
ysis, computer  programs  were  developed  to  provide  for  computerized 
operation  of  the  various  models. 

Ultimately,  deterioration  and  reliability  models  for  the  rutting 
mode  were  developed  for  unsurfaced,  gravel- surfaced,  and  two-  and  three- 
layer  flexible  pavements,  respectively.  These  models  were  then  combined 
to  provide  for  deterioration  and  reliability  analysis  as  well  as  for  dif- 
ferential damage  analysis  where  mixed  traffic  occurs.  This  system  is 
termed  the  Differential  Analysis  System  (DAS). 

Definition  of  Terms 

For  clarity,  certain  terms  pertinent  to  this  document  are  listed  and 
defined  as  follows: 

1.  Pavement . — A horizontal  structure  intended  to  protect  a subgrade 
from  the  loading  effects  of  wheeled  vehicles. 


asphalt  concrete  surface  course,  base  course,  and  subbase  course  above 
the  subgrade. 


3.  Two-layer  flexible  pavement. — A pavement  consisting  of  an  as- 
phaltic concrete  surface  course  and  a base  course  above  the  subgrade. 

Gravel-surfaced  facility. — Any  facility  intended  for  use  of 
wheeled  vehicles  and  where  a gravel  course  serves  as  the  pavement 
structure. 

5.  Unsurfaced  facility. — A facility  intended  for  use  of  wheeled 
vehicles  and  where  no  pavement  structure  exists  above  the  in  situ 
material. 


6.  Equivalent  single-wheel  load  (ESWL). — That  load  on  a single 
wheel  that  produces  the  same  effect  (usually  measured  in  terms  of  ver- 
tical deflection)  beneath  the  wheels  as  a group  of  wheels  with  the  same 
single-wheel  contact  area. 

7.  Operation  or  repetition. — One  pass  of  one  vehicle  over  a section 
of  pavement. 

8.  Coverage. — One  pass  of  a wheel  over  every  point  in  a trafficable 


area. 


9*  Serviceability. — The  capability  of  a facility  to  perform  the 
intended  functions. 

10.  Deterioration. — Any  departure  from  the  as-constructed  condition 
of  a facility  that  results  in  a reduction  in  serviceability. 

11.  Deterioration  mode. — The  nature  of  deterioration. 

12.  Structural  mode. — A mode  of  deterioration  that  is  in  terms  of 
the  structural  properties  or  capabilities  of  a facility. 
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13.  Functional  mode — A mode  of  deterioration  that  is  in  terms  of 


a reduction  in  serviceability. 

1^.  Life-cycle  management Quantative  optimization  of  design, 

maintenance,  and  repair,  with  respect  to  serviceability  of  a facility 
from  conception  to  the  end  of  its  life. 
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DESIGN  PROCEDURES  AND  PREVIOUS  RESEARCH 

Design  Procedures 

Several  design  procedures  currently  exist  for  the  determination  of 
thickness  and  strength  requirements  for  protection  of  subgrades.  Some 
designs  currently  in  use  are  the  Texas,  CBR  (CE),  Group  Index,  Califor- 
nia, FAA,  and  Asphalt  Institute  methods.  These  designs  are  empirically 
based,  theory  based,  or  in  some  cases  a combination  of  both.  The  advent 
of  computers  brought  about  more  extensive  utilization  of  elastic  theory 
in  developing  theoretical  design  procedures.  Among  these  are  the  Shell, 
Chevron,  Asphalt  Institute,  and  the  recently  developed  CE  design 
procedures. 

These  procedures  have  served  the  respective  agencies  well  in  the 
functions  intended.  However,  as  a general  rule  they  have  been  applicable 
only  to  preconstruction  design  of  pavements.  Some  of  the  agencies  have 
modified  their  procedures  for  use  in  pavement  condition  surveys  and  for 
overlay  design  tools. 

The  past  decade  has  seen  record  pavement  construction  of  all  types. 
Many  of  these  pavements  either  are  approaching  or  have  already  exceeded 
their  respective  design  lives  and,  therefore,  lie  in  some  state  of  de- 
terioration and  need  of  repair.  The  design  procedures  previously  men- 
tioned have  usually  been  found  inadequate  as  tools  for  quantification 
of  deterioration  and  subsequent  repair  needs.  Therefore  an  era  has 
arrived  where  maintenance  and  rehabilitation  are  at  the  forefront  and 
where  quantitative  analysis  procedures  are  either  nonexistent  or 
lacking  in  adequacy  or  validation. 

The  CE  design  procedures  for  rigid  (l6),  flexible  (13),  and 
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gravel-surfaced  pavements,  as  well  as  unsurfaced  facilities  (12),  fall 
into  such  a category  and  serve  as  examples  of  design  procedures  that  do 
not  provide  for  analysis  of  deterioration  or  the  reliability  of  a 
pavement . 

Deterioration  and  Reliability  Investigations 

The  concepts  of  statistics  and  probability  are  certainly  not  new  to 
the  field  of  engineering  as  a whole.  However,  the  most  significant 
inroads  to  the  utilization  of  these  concepts  have  been  made  in  recent 
years.  The  First  International  Conference  on  Applications  of  Statistics 
and  Probability  to  Structural  and  Soils  Engineering  was  held  in  Hong  Kong 
in  1971  (22).  Material  presented  at  that  conference  represents  signifi- 
cant beginnings  in  the  overall  application  of  these  tools  in  the  field 
of  civil  engineering. 

The  utilization  of  probability  and  statistics  to  address  deteriora- 
tion and  reliability  of  pavement  life-cycle  management  is  newer  still. 
Significant  contribution  was  made  in  197^  by  Lu,  Lytton,  and  Moore  (20) 
for  the  Texas  Transportation  Institute  in  cooperation  with  the  Federal 
Highway  Administration.  They  utilized  data  collected  from  pavement  test 
sections  in  Texas  to  predict  serviceability  loss  in  flexible  pavements. 
The  concept  of  probabilistic  design  used  was  formulated  by  Darter  and 
Hudson  and  applied  to  flexible  pavement  design  systems  in  1973  (9).  Lu, 
Lytton,  and  Moore  developed  a two-step  constrained  select  regression  pro- 
cedure to  examine  the  effect  of  each  variable  on  pavement  serviceability 
loss.  They  also  used  stochastic  reliability  concepts  to  evaluate  ex- 
pected value  and  variance  of  serviceability  loss.  In  1975.  Hudson  et  al. 
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(15,17) , in  a contract  study  for  the  WES,  wherein  the  state  of  the  art 
in  predicting  pavement  reliability  was  reviewed,  recommended  that  re- 
search should  continue  in  the  area  of  pavement  reliability.  The  current 
status  of  reliability  assessment  was  also  given  in  1975  by  Barker  and 
Brabston  {h) . Although  the  new  design  procedure  for  flexible  airport 
pavements  (4)  is  innovative  and  provides  greater  capabilities  in  design, 
they  state  that  pavement  deterioration  is  a continuous  function  but  is 
treated  as  discontinuous  by  criteria  that  label  pavements  as  "failed"  or 
"unfailed."  They  further  state  that  this  is  not  the  case  but  that  un- 
fortunately methodology  still  does  not  exist  to  predict  deterioration 
realistically. 

The  net  result  of  these  investigations  is  that  not  only  is  sto- 
chastic reliability  and  deterioration  a viable  approach  to  life-cycle 
management  in  view  of  material  variability  and  other  uncertainties,  but 
the  state  of  the  art  exists  for  the  application  of  these  concepts  to 
various  design  procedures,  such  as  the  CE  design  procedure  as  illustrated 
in  this  document. 

The  Pavement  Deterioration  Program 

The  need  to  more  effectively  construct  new  pavements  and  the  neces- 
sity to  maintain  many  existing  pavements  has  been  a concern  of  various 
Federal  agencies.  A pavement  deterioration  and  reliability  analysis 
program  was  instituted  at  the  WES  in  1974  and  sponsored  by  various 
agencies  for  the  purpose  of  developing  methodology  for  the  effective 
life-cycle  management  of  pavements.  As  the  initial  achievement, 
methodology  has  been  developed  to  analyze  deterioration  and  to  assess 
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reliability  in  terms  of  rutting.  These  concepts  have  been  developed  to 
be  applicable  to  the  original  CE  design  procedure  for  flexible  pavements 
and  to  utilize  the  parameters  of  that  procedure.  Rutting  was  chosen  as 
the  initial  mode  of  deterioration  in  which  to  test  the  hypothesis  that 
such  achievements  could  be  made  due  to  availability  of  data  and  since 
rutting  has  historically  constituted  failure  criteria  for  CE  flexible 
pavements . 

This  document,  in  the  succeeding  parts,  will  describe  the  test  pro- 
gram as  well  as  the  deterioration  and  reliability  assessment  methodology. 

The  overall  research  and  data  collection  program  is  aimed  at  deteri- 
oration and  reliability  assessment  for  rigid,  flexible,  and  all  other 
types  of  pavements  and  includes  all  modes  of  deterioration.  However, 
this  research  is  intended  to  set  the  technological  framework  for  over- 
all analysis  by  providing  the  procedures  for  dealing  with  all  types  of 
deterioration  on  all  pavements.  When  these  procedures  are  eventually 
applied  to  all  deterioration  modes  and  all  pavements,  a complete  life- 
cycle  management  system  can  be  developed. 
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DATA  COLLECTION  AND  CURRENT  STUDIES 

Earlier  Tests  and  Data  Collection 

The  CE  design  method  has  required  revalidation  and  revision  since 
its  adoption  due  to  the  ever-changing  nature  of  traffic.  Airfield  design 
criteria  have  seen  the  greatest  change  due  to  the  increase  in  aircraft 
weights,  wheel  loads,  tire  pressures,  and  number  of  operations.  Road 
and  highway  design  criteria  have  also  changed  considerably  over  the  years 
due  to  changes  in  vehicle  characteristics  and  modes  and  quantities  of 
operations.  To  stay  abreast  of  these  changes  and  to  provide  the  most 
applicable  criteria  possible,  it  has  been  necessary  for  the  CE  to  con- 
duct numerous  prototype  tests.  These  prototype  tests  have  classically 
proved  the  best  approach  until  recently  when  such  studies  have  become 
cost  prohibitive.  The  net  result  of  such  an  extended  series  of  tests 
has  been  the  accumulation  of  myriad  prototype  pavement  performance  data. 

The  data  are  necessarily  in  terms  usable  in  the  early  CE  design  methods , 
namely  rut  depths,  thickness  of  layers,  strengths  of  layers  and  subgrades 
in  terms  of  CBR,  and  vehicle  characteristics,  which  include  number  and 
configuration  of  wheels,  tire  pressures,  and  wheel  loads.  Although  the 
data  were  accumulated  under  closely  controlled  conditions,  the  re- 
searchers were  attempting  to  determine  end-point  fa:  Lure.  The  design 
method  was  structured  to  determine  thicknesses  required  to  prevent  sub- 
grade deformation  as  a result  of  loading,  and  the  failure  criterion  was 
largely  in  terms  of  maximum  allowable  rut  depths.  Therefore,  in  the 
process  of  repetitive  loading  and  intermittent,  frequent  data  collection, 
large  quantities  of  data  were  collected  that  characterize  a change  in 
rut  depth  as  a function  of  traffic. 
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Deterioration  Data  Search 


A data  search  and  analysis  was  conducted  in  the  early  stages  of  the 
deterioration  program  to  determine  whether  suitable  data  existed.  The 
study  consisted  of  a two  man-year  effort  to  review  all  existing  data  at 
the  WES  and  to  screen  the  data  for  overall  applicability  prior  to  initi- 
ating the  development  research.  The  results  were  that  a large  quantity 
of  high-quality  rutting  data  was  available  that  applied  to  two-  and 
three-layer  flexible  pavements  as  well  as  gravel- surfaced  and  unsurfaced 
facilities.  The  data  were  initially  termed  suitable  if  the  variables 
were  within  appropriate  ranges  and  were  all  included  and  deliberately 
recorded.  Further  analysis  of  distribution  and  range  of  data  is  pre- 
sented later  in  this  document. 

The  literature  reviewed  and  utilized  as  data  sources  is  tabulated 
for  reference  uses  in  Appendix  II.  The  90  references,  which  are  not 
referred  to  individually,  are  considered  to  comprise  results  of  essen- 
tially all  key  tests  at  the  WES  in  recent  years. 

A tabulation  of  the  data  ultimately  used  in  this  study  is  given  in 
Appendix  III.  These  data  represent  the  final  data  selected  for  develop- 
ment of  deterioration  and  reliability  models  and  are  by  no  means  the 
total  data  available.  Criteria  for  rejection  of  certain  data  were  gen- 
erally based  on  range  and  reaction  of  test  pavements.  Rejection  cri- 
teria will  be  discussed  along  with  model  development. 

A prototype  pavement  test  section  typical  of  those  upon  which 
traffic  tests  were  conducted  and  data  collected  is  shown  in  Fig.  1.  Such 
test  areas  are  constructed  under  rigid  controls  to  a specified  design 
requirement.  Traffic  is  normally  applied  in  a specified  pattern  to 
simulate  actual  distribution.  Specially  designed  trafficking  vehicles 
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having  load  wheels  that  exactly  represent  the  desired  tire  sizes,  pres- 
sures, loads,  and  configurations,  as  shown  in  Fig.  2,  are  employed  to 
place  the  desired  traffic  upon  the  test  area. 

Various  materials  tests,  surveys,  and  surface  measurement  are  made 
prior  to,  during,  and  after  traffic  to  monitor  all  physical  conditions 
and  all  changes.  Fig.  3 gives  a typical  rut  depth  measurement,  while 
Fig.  4 shows  a strength  (CBR)  measurement  in  progress  on  one  of  the 
layers. 

Although  Figs.  1-U  illustrate  typical  prototype  test  sections  built 
and  tested  to  generate  much  of  the  data,  a significant  quantity  of  data 
was  accumulated  on  actual  roads  and  airfields  that  were  subjected  to  de- 
sign traffic  either  on  an  accelerated  basis  or  in  the  normal  mode  of 
operation. 

The  availability  of  such  a quantity  of  deterioration  data  in  terms 
of  rutting  on  CE-designed  pavements  is  the  reason  for  selecting  rutting 
and  CE  design  parameters  for  use  in  this  basic  research  on  deterioration 
and  reliability.  As  data  are  made  available  under  the  test  program 
described  in  succeeding  paragraphs,  the  developments  reflected  in  this 
document  will  be  expanded  to  other  deterioration  modes  and  to  other  de- 
sign procedures  according  to  the  needs  of  other  governmental  agencies. 

Current  Data  Collection  Programs 

The  CE-FS  agreement  (10)  called  for  the  pilot  studies  to  develop 
methodology  as  reported  herein  and  for  a large-scale  field  testing  pro- 
gram. The  purpose  of  the  program  was  to  accumulate  actual  deterioration 
data  pertinent  to  FS  and  CE  roads  for  use  in  validation  of  this  procedure 
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for  employment  by  the  respective  agencies.  The  modes  of  deterioration 
considered  paramount  for  this  study  were  rutting,  roughness,  slipperi- 
ness, cracking,  and  surface  loss  on  aggregate- surfaced  roads.  These 
field  tests  have  been  in  progress  for  approximately  two  years.  Initial 
liaison  was  established  in  most  cases  by  FS  personnel  who  also  assisted 
in  establishing  test  sites  at  several  locations  throughout  the  United 
States.  Currently,  approximately  hO  test  sites  exist  that  are  respec- 
tively being  monitored  for  deterioration  in  the  various  modes.  It  is 
anticipated  that  the  program  will  continue  for  approximately  four  more 
years.  During  the  latter  stages,  data  collected  will  be  applied  to  the 
system  developed  herein. 

Test  sites  have  been  selected  at  the  various  regional  locations  in 
areas  having  suitable  design  and  traffic  features  to  provide  a deterio- 
ration environment.  Test  sections  are  established  at  these  test  sites 
whereon  surface  conditions  and  pavement  layer  strengths  and  thicknesses 
are  monitored  periodically.  Fig.  5 gives  a layout  used  at  most  test 
sites.  The  layout  shows  locations  of  test  pits  as  well  as  locations  for 
profiling  for  roughness  and  cross-section  and  rut  depth  measurements. 
These  tests  are  conducted  in  a conventional  manner  and  in  sufficient  de- 
tail to  depict  any  deterioration.  In  addition  to  these  tests,  roughness 
is  monitored  using  a Mays  ride  meter.  Skid  resistance  is  measured  in 
terms  of  energy  required  to  produce  slip  or  loss  of  traction.  The  slip- 
energy  device  and  recorder  were  especially  designed  for  this  test.  The 
data  collected  from  this  test  are  also  expected  to  be  applicable  to  slip- 
energy  studies  being  conducted  by  the  FS  (9)> 

Fig.  6 illustrates  rut  depth  measurement  on  a one-lane  gravel- 
surfaced road.  The  device  used  is  a standard  10- ft  straightedge,  and 
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-TYPICAL  TEST  SITE  LAYOUT 


the  rut  depth  is  considered  to  be  the  maximum  deviation.  Fig.  7 shows 
survey  personnel  conducting  profile  and  cross-section  measurements,  both 
of  which  are  normally  taken  at  1-ft  intervals  to  provide  for  the  maximum 
practical  definition  of  the  surface  configuration.  Fig.  8 shows  small 
aperture  testing  (SAT)  of  layer  strengths.  The  SAT  procedure  (11)  pro- 
vides a 6-in.-diam  access  hole  through  which  layer  strengths  CBR  values 
as  well  as  thicknesses  in  inches  can  be  determined.  Additionally,  sam- 
ples can  be  retained  for  moisture  content  determination.  Moisture 
content  determination  is  frequently  made  using  nuclear  testing  devices 
(26)  as  is  illustrated  in  Fig.  9.  Such  procedures,  when  used  in  con- 
junction with  SAT,  provide  for  more  rapid  and  economical  monitoring  of 
changes  in  pavement  conditions. 

In  addition  to  these  tests,  climatological  data  are  collected  for 
the  area  on  a continuous  basis. 

Monitoring  to  determine  type  and  quantity  of  traffic  that  brings 
about  deterioration  has  been  of  primary  concern  and  interest  throughout 
the  program.  Several  procedures  and  items  of  equipment  have  been 
utilized  with  varying  degrees  of  success.  However,  the  greatest  success 
to  date  has  been  by  use  of  an  inductive  loop  counter  synchronized  with 
a 35-mm  movie  camera.  This  combination  provides  for  not  only  a traffic 
count  but  a sequence  of  photographic  frames  that  depicts  each  vehicle 
crossing  the  loop,  including  those  conducting  maintenance.  Fig.  10 
shows  the  manner  of  installing  the  induction  loop,  while  Figs.  11  and  12 
show  the  hidden  counter  and  camera,  respectively.  Fig.  13  gives  an  ex- 
ample of  heavy  traffic  that  frequents  many  roads,  especially  in  timber 
sale  operations  in  national  forests.  Figs.  1^  and  15,  show  maintenance 
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■RUT  DEPTH  MEASUREMENT  AT  A SITE  ON  A GRAVEL- 
SURFACED ONE-LANE  FOREST  ROAD 


FIG.  7.— SURFACE  CHARACTERIZATION  WITH  CROSS-SECTION 
AND  PROFILE  MEASUREMENTS  AT  A TEST  SITE 
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FIG.  10.— MANNER  OF  INSTALLATION  OF  INDUCTION  LOOP 
FOR  TRAFFIC  COUNTING 


FIG.  11.— VIEW  OF  HIDDEN  TRAFFIC  COUNTER 
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VIEW  B 


VIEW  A 


FIG.  12.— VIEW  OF  HIDDEN  CAMERA 


FIG.  15 WATERING  IS  NECESSARY  TO  MINIMIZE 

DUST  PROBLEMS  ON  FEEDER  ROADS 


operations  on  an  unsurfaced  logging  road.  Use  of  the  induction  loop 


counter  and  camera  will  not  only  aid  in  counting  and  describing  traffic 
but  also  serve  as  a permanent  maintenance  record. 

Tests  similar  to  those  described  above  are  being  conducted  at  mili- 
tary posts  and  in  private  industry  forests  in  the  South  in  cooperation 
with  the  respective  agencies  in  an  attempt  to  further  expand  the  data 
base.  The  combination  of  tests  and  traffic  monitoring,  when  successfully 
pursued  over  a period  of  approximately  six  years,  is  expected  to  provide 
abundant  data  for  deterioration  and  reliability  system  validation. 

Future  Data  Collection 

Plans  have  tentatively  been  formulated  to  expand  the  testing 
procedure  to  include  aggregate  surface  loss  tests.  The  surface  loss  mode 
of  deterioration  is  critical  due  to  the  high  cost  of  replacement  of 
surface  aggregate  lost  each  year  as  a result  of  both  traffic  and  mainte- 
nance operations.  The  initial  stages  of  surface  loss  studies  would  con- 
sist of  development  of  test  procedures  as  well  as  preliminary  determina- 
tion of  pertinent  variables.  Although  aggregate  loss  is  not  held  in 
high  regard  universally  as  a deterioration  mode,  it  carries  significant 
impact  with  the  FS.  It  has  been  determined  that  annual  cost  of  replace- 
ment of  aggregate  on  gravel- surfaced  roads  is  a major  cost  item  for  that 
agency.  Inroads  have  been  made  into  the  surface  loss  problem  by  way 
of  slip-energy  and  associated  tire  wear  concepts  (9).  However,  it  is 
hypothesized  that  another  valid  approach  is  the  deterioration  and  reli- 
ability concept  set  forth  in  this  document.  The  cooperating  agencies 
therefore  hope  to  achieve  results  through  at  least  one  of  the  approaches. 
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The  thrust  of  these  paragraphs  has  been  to  set  forth  the  overall 
philosophy  of  the  program.  Namely,  in  order  for  the  overall  work  to  be 
accomplished  in  the  foreseeable  future,  existing  data  are  used  to  explore 
the  hypothesis  that  deterioration  and  reliability  concepts  can  be  em- 
ployed to  effect  life-cycle  management  while  a full-blown  data  collection 
program  for  validation  is  in  progress.  The  succeeding  portion  of  this 
document  reports  findings  as  to  the  development  of  pilot  procedures  for 
predicting  the  deterioration  and  reliability  of  pavements. 
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DEVELOPMENT  OF  DETERIORATION  MODELS 


Initial  Data  Analysis 

The  primary  objective  in  examining  the  hypothesis  that  the  rate  of 
deterioration  of  a pavement  structure  can  be  quantified  was  the  estab- 
lishment of  the  rate  of  rut  depth  (RD)  change  as  a function  of  the  inde- 
pendent variables.  In  the  case  of  this  work,  where  up  to  eight  inde- 
pendent variables  were  involved,  it  was  recognized  early  that  the 
method  of  data  analysis  would  be  critical  in  terms  of  time,  cost,  and 
overall  results  of  the  analysis.  The  initial  analysis  consisted  of 
utilizing  a conventional  one-step  regression  for  the  entire  set  of  vari- 
ables involved.  The  procedure  quickly  proved  to  be  rather  time-consuming 
and  produced  rather  poor  correlation  in  cases  where  it  was  intuitively 
obvious  that  better  results  could  be  attained.  With  this  experience  in 
hand,  other  procedures  were  considered  for  analysis,  including  impos'ng 
forms  and  coefficients  onto  the  variables  and  performing  regressions  on 
these  "new"  variables.  Again,  the  results  were  not  satisfactory. 

The  initial  data  analysis,  although  largely  unsuccessful,  provided 
at  least  two  results  that  would  prove  valuable  in  future  efforts. 

First,  it  was  determined  that  some  reasonable  correlations  did  exist 
between  the  rate  of  rutting  and  the  nature  of  traffic  on  given  pavement 
structures.  Second,  it  became  possible  to  make  the  judgment  that  im- 
posing what  appeared  to  be  a suitable  form  on  the  variables  in  some  in- 
stances could  effectively  provide  for  a better  overall  correlation  and  a 
more  direct  route  in  attaining  that  correlation. 
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Regression  Analysis  Procedures 


The  procedure  of  analysis  ultimately  selected  for  use  is  an  orderly 
method  of  developing  mathematical  relations  from  sets  of  data  using  mul- 
tiple regression  analysis.  This  procedure  was  developed  by  Lu,  Lytton, 
and  Moore  (20),  although  the  basic  method  has  probably  been  employed  by 
others,  due  to  its  direct  approach  to  the  analysis  of  data.  The  pro- 
cedure was  first  described  in  reference  20,  although  it  was  first  used 
by  Lytton  and  Castleberry  (6)  in  a study  on  damage  to  houses  located  on 
expansive  clays. 

As  originally  developed,  the  procedure  used  the  SELECT  regression 
program,  which  had  been  developed  at  Texas  A&M  University  by  researchers 
with  the  Institute  of  Statistics,  namely  R.  R.  Hocking,  R.  N.  Leslie, 

L.  R.  La  Motte,  and  D.  A.  Debuse.  Their  development  of  the  SELECT  re- 
gression method  and  computer  program  is  recorded  in  references  8,  lU, 
and  19*  The  SELECT  procedure  is  desirable  because  of  its  several  good 
features  and  was  utilized  in  the  earlier  stages  of  development  of  equa- 
tions shown  herein.  However,  the  final  regressions  performed  herein 
were  conducted  using  a WES  library  program  called  STAT21.  The  program 
is  similar  to  SELECT  and  was  confirmed  to  produce  nearly  identical  re- 
sults. Its  availability  at  WES  dictated  its  use  for  this  study. 

The  procedure  utilized  is  termed  a two-step  constrained  regression 
procedure.  The  first  step  of  this  method  is  essentially  a selection  re- 
gression procedure  using  a multiplicative  model  to  obtain  the  approximate 
exponents  of  the  independent  variables.  The  second  step  determines  the 
coefficients  of  linear  combinations  of  the  products.  The  final  model  is 
selected  based  upon  four  factors: 
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1.  As  simple  an  expression  as  can  reasonably  describe  what  is 
actually  taking  place  in  terms  of  the  dependent  variable  and  reflecting 
the  effect  of  the  independent  variables. 

2.  High  multiple  correlation  of  the  model. 

3.  Small  prediction  error. 

U.  Satisfaction  of  physical  constraints. 

The  procedure  used  in  this  study  represents  an  alteration  of  the 
original  procedure  utilized  in  program  SELECT.  The  departure  was  pos- 
sible due  to  the  fact  that  normally  there  were  several  measured  values  of 
rut  depth  that  increased  according  to  the  number  of  load  repetitions  for 
each  series  of  data.  Since  the  data  were  composed  of  rut  depth-load  his- 
tories, the  regression  itself  was  broken  into  several  stages  as  described 
in  the  following  paragraphs. 

The  variables  involved  in  the  equations  are  tabulated  as  follows: 

RD  = rut  depth,  in. 

P = equivalent  single-wheel  load,  lb 

t = tire  pressure,  psi 
R = number  of  load  repetitions 

t^  5 layer  thickness,  in. 

= layer  strength  in  CBR 

The  total  number  of  variables  changes  from  one  pavement  type  to 
another  depending  upon  the  number  of  layers  that  make  up  the  payment. 

The  procedure  described  below  is  in  general  terms  and  applicable  to  all 
of  the  pavement  types. 

The  first  step  in  equation  development  is  determining  the  rut  depth 
as  a function  of  the  traffic  variables,  which  are  equivalent  single-wheel 


29 


load,  tire  pressure,  and  number  of  repetitions.  The  equations  thus  de- 
rived will  take  the  form: 


RD  = ^(P.tp.R)  + 

where  SE^  is  the  regression  error  due  to  lack  of  fit  and  to  variables 
that  have  not  yet  been  included. 

The  next  step  was  a regression  upon  a "normalized"  rut  depth,  NRD, 
given  by  the  following  equation: 


NRD 


RD 

(P,tp,R) 


f2<Vi>  * SE2 


This  is  termed  "normalized"  because  any  effects  of  loading  are  accounted 
for  in  the  new  dependent  variable.  The  remaining  variables  (pavement 
structure)  were  used  to  show  how  their  properties  influenced  rut  depth. 
The  new  regression  error,  SE0  , is  affected  by  three  factors,  which  are 
lack  of  fit  of  the  new  equation,  stochastic  variation  of  the  variables 
and  climatic  variations,  and  other  factors  not  explained  in  the 
variables. 

The  third  step  consisted  of  multiplying  f by  f0  and  running  the 
regression  a third  time  to  correct  the  nonlinearities  that  were  intro- 
duced in  the  first  two  steps.  The  form  of  this  third  equation  becomes 

RD  = f + SE.  = at^fjh  + SE3 


This  step  is  the  final  regression  step  and  generally  results  in  a higher 
correlation  coefficient,  r , and  lower  error,  SE,  , than  all  the 
previous  models. 
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The  last  step  is  to  determine  the  size  of  the  error,  SE^  , which  is 
accomplished  by  getting  the  sum  of  square  errors,  a2  , between  the  ob- 
served and  mean  rut  depths.  This  step  is  expressed  as 

2 A (RD.  - RD. )2 
at = 

i=l 

where 

RD  = the  mean  rut  depth 
RD  = the  observed  value  of  rut  depth 
This  sum  of  square  errors  is  composed  of  the  sum  of  the  variances  due  to 
the  regression  equation  used  and  the  lack  of  fit  of  the  equation 


where 


Jtof 

The  error 


by 


= variance  due  to  regression 
= variance  due  to  lack  of  fit 

SE^  is  measured  by  the  lack-of-fit  variance  and  is  determined 


2 

atof 


^ (RD.  - RD.)2 
i=l 


n - 1 


where  RDi  is  the  value  of  the  rut  depth  predicted  by  the  regression 


equation.  The  variance  due  to  regression  is  determined  by 

2 


n (RD.  - RD)2 


2 - kd, 

3r  Zj  n - 1 


i=l 


or  the  regression  variance  may  be  approximated  by  a method  based  upon  the 
Taylor's  series  expansion  of  the  predictive  function. 
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Development  of  Regression  Models 


The  data  shown  in  Appendix  II  were  used  to  develop  models  for 
rutting  in  the  manner  described  in  the  preceding  paragraphs.  Experience 
and  familiarity  in  the  area  of  CE  pavement  design  procedures  provided 
some  insight  into  the  probable  behavior  of  the  variables;  therefore,  some 
transformations  were  attempted  in  the  case  of  many  of  the  variables.  The 
thickness  equations  currently  in  use  by  the  CE  (13)  were  studied  to  de- 
termine the  most  likely  form  that  some  of  the  variables  might  assume. 
Although  the  procedure  essentially  provides  for  determining  thickness 
required  to  protect  subgrades,  rutting  was  normally  the  limiting  cri- 
terion. Therefore,  the  design  equations  can  be  considered  as  having  the 
ability  to  predict  one  certain  value  of  a rut  depth.  If  an  inversion  of 
variables  is  envisioned,  one  can  see  that  even  in  those  equations, 
rutting  could  be  considered  "dependent"  and  thickness  could  be  a con- 
trolled variable.  If  this  were  the  case,  repetitions  could  be  varied  to 
determine  thickness  required  to  produce  the  limiting  rut  depth.  Thus, 
certain  forms  were  imposed  upon  the  variables  in  an  attempt  to  intui- 
tively achieve  better  representation  of  the  data  in  the  regression 
analysis. 

First  among  these  considerations  was  the  indication  from  sample 
data  that  the  rut  depth  increased  in  some  cases  with  the  square  root  of 
tire  pressure  and  in  others  with  the  common  logarithm  of  tire  pressure. 
Therefore,  these  two  forms  were  attempted  and  can  be  seen  in  the 
equations. 

Additionally,  thickness  of  some  layers  seemed  to  have  various  types 
of  correlation  with  rut  depth.  Some  examples  of  the  forms  for  thickness 
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and  tire  pressure  that  were  intuitively  developed  and  employed  in  the 
final  equation  are  as  follows: 


RD  = f(t  )1//2 

P 

RD  = f[log(t  )]  ... 

RD  = f[l  + log(t1)  ] 

RD  = fCa^  + bt2)  . 


where  a and  b were  constants  usually  having  values  of  1.00  and  1.25, 
respectively,  which  are  weighted  coefficients  intuitively  selected. 

The  following  four  parts  of  this  document  give  a description  of  the 
development  of  the  models  for  the  four  categories  of  pavements  where 
applicable  data  existed. 

Unsurfaced  Facility  Model 

The  data  available  for  development  of  rutting  models  fell  into  two 
categories.  The  first  category  is  one  where  the  surface  CBR  (C^)  is 
usually  equal  to  or  lower  in  value  than  the  subgrade  CBR  (C2).  The  sec- 
ond category  is  when  is  greater  than  C2  and  will  be  discussed  sub- 

sequently. These  data  typify  many  facilities  where  in  situ  materials  are 
used  and  the  facility  is  constructed  merely  by  grading  a smooth  surface. 
The  CBR  in  the  upper  portion  is  usually  lower  due  to  moisture  or  organic 
materials,  while  the  lower  (subgrade)  portion  is  more  stable  due  to 
lesser  moisture  fluctuation  and  other  disturbances. 

The  data  in  this  category  were  inspected,  and  a total  of  1^2  data 
points  were  retained  for  analysis.  The  variables  influencing  the  rut 
depth  are  shown  in  Table  1. 
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TABLE  1.— VARIABLES  FOR  UNSURFACED  AND  GRAVEL-SURFACED 
FACILITY  RUTTING  MODELS 


* C^  < Cg  for  unsurfaced  facilities. 

C^  > C0  for  gravel-surfaced  facilities. 

The  data  used  in  development  of  the  unsurfaced  model  are  given  in 
Appendix  III.  Numerous  iterations  were  performed  using  the  regression 
analysis  techniques  previously  described.  As  models  were  developed,  they 
were  rejected  on  one  or  more  of  the  four  criteria  previously  discussed. 
The  usual  basis  for  rejection  was  low  correlation  of  the  data.  However, 
several  models  were  rejected  due  to  improper  behavior  of  variables  or 
simply  to  unseemly  appearance  of  the  model  when  compared  with  forms  found 
in  the  earlier  pavement  thickness  design  equations. 

Of  the  lb2  unsurfaced  data  points  where  the  surface  CBR  was  lower 
than  the  subgrade  CBR,  8 points  that  had  rut  depths  9 in.  or  higher  were 
eliminated  since  it  was  considered  unreasonable  to  predict  such  high 
values.  For  the  13^  remaining  points,  the  equation  considered  best  was 


3U 


0.00609  tt- 
L2 


p 0.U336  1Q46L  .1.0670  0. 5226 

— t (log  t)  R 


,2.0267 


where  the  standard  error  of  estimate  (SE)  and  the  coefficient  of  corre- 
lation (r)  equal  0.530  and  0.9137,  respectively.  This  model  shows  an  in- 
crease in  RD  with  an  increase  in  P , t , t , and  R and  a decrease 

P 

in  RD  with  an  increase  in  and  Cp  . It  should  be  noted  that  in 

the  case  of  the  gravel-surfaced  facility  model  where  the  top  layer  was 
stronger,  the  RD  shows  a decrease  with  higher  thicknesses,  but  in  the 
data  points  where  the  top  layer  is  weaker,  the  RD  shows  an  increase 
with  higher  thicknesses.  For  this  equation,  the  exponent  of  is  low 

compared  with  the  exponent  of  . However,  this  equation  was  con- 
sidered the  best  by  comparison.  Using  the  residuals  of  this  equation, 

20  "bad"  data  points  were  deleted  leaving  11*4  points. 

Using  the  remaining  11 14  data  points,  the  equation  considered  best 
was  derived  using  the  step  regression  as  follows: 

1.  Step  1. — Solve  for  RD  as  function  of  four  variables: 


RD  = 0.0025 


Oo.  14976  0. 8628  0. 5070 

run 

E „ 

„1.9976 


where  SE  = O.I4O8  and  r = 0.9373 

2.  Step  2. — Divide  RD  by  the  expression  on  the  right,  and  solve 
for  the  other  two  variables: 


p0. U976  t0.8628  ^0. 5070 

^1.9976 


= 0.108 


(log  t 


,0.14337 


,1. 213c 
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3.  Step  3. — Multiply  expression  on  right  by  denominator  on  left: 


RD  = 0.110 


p0. 1+976  ^ 0.8628  R0.  5070  (lQg  t)O.U337 

P -1.9976  -1.2138 
L1 


0.9898 


where  SE  = 0.399  and  r = O.9I+O3 


1+.  Step  k. — Multiply  exponent  of  each  variable  by  O.9898,  and  the 
final  equation  becomes 


RD  = 0.110 


p0. 1+925  ,0.851+8  R0.  5018  (logt)0.U293 


,1-9773  -1.2015 
'1  L2 


where  SE  = 0.399  and  r = O.9I+03 


With  this  step  the  outside  radical  can  be  eliminated. 

This  is  the  model  finally  selected  as  a result  of  regression 
analysis.  Other  equations  and  partial  step  equations  that  were  tried 
and  eliminated,  along  with  the  reasons  for  eliminating  them,  were  as 
follows : 


RD  = 17.773 


-0.1+712  1.0005  ^0.131+9  D0. 1+610 

r t t n 

E 

1. 9899  2. 965I 

L1  ^2 


where  SE  = 0.1*01  and  r = 0.9391* 


exponent  of  t too  low;  exponent  of  Cp  too  high 
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r 


RD  = 0.0214 


0.5042 


t0.8895  R0. 4889 
.0.1941  .2.011 


where  SE  = 0.406  and  r = 0.9380 
not  in  numerator 


RD  = 0.00254 


p \0.1038 


1.6943 


rO.5545 

to. 0765 


where  SE  = 0.557  and  r = 0.8797 

p 

t^  not  in  numerator;  exponent  of  — too  low 


RD  = 9,383,500 


1.7^52 


p0.l490  ^0. 9184  R0.44l7 

.8.029 

^2 


where  SE  = 0.492  and  r = 0.9076 

coefficient  too  high;  exponent  of  Cp  too  high; 

exponent  of  P too  low 


p0.5053  t0.8692  r0.4966 

RD  = 0.00296  — ■■  ■ „ P 7r-s-A‘b 

c^.ono  (log  t) 0.6198 

where  SE  = 0.409  and  r = 0.9371 

log  t not  in  numerator 
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RD  = 26.302 


pO.UTOO  ,1.0025  (lQg  t)0.4220  RQl*6l3 


where  SE  = 0.1*01  and  r = 0.9395 


exponent  of  too  high 


Gravel-Surfaced  Facility  Model 

The  gravel-surfaced  facility  data  include  a total  of  299  data 
points.  The  data  typify  a facility  where  a gravel  or  similar  surface 
generally  referred  to  as  "gravel"  is  placed  upon  the  existing  subgrade 
to  provide  for  protection  of  that  subgrade.  Therefore,  these  data  are 
characterized  by  gravel-surfacing  CBR  values  (C^)  higher  than  the  CBR 
values  of  the  subgrade  (C^).  The  variables  shown  in  Table  1 are  appli- 
cable to  this  model  also.  However,  it  should  be  noted  that  C repre- 
sents a surfacing  of  some  type  of  gravel  with  a finite  thickness  t . 

Numerous  regressions  were  performed  on  the  original  299  data  points, 
and  the  final  equation  was  derived  in  steps  as  follows: 


1.  Step  1. — Rut  depth  was  determined  as  a function  of  three 


variables : 


/ \0.5U39  R0. 161*7 

RD  = 0.0123  £- 

\C1 

where  SE  = 0.6ll  and  r = 0.6809 


2.  Step  2. — These  three  variables  were  transposed,  and  regression 
was  performed  on  the  remaining  three  variables: 


I 


RD 

0.51+39 


R0. 16U7 


0.0103 


t0.  53l*3 

P 

,0.6911*  0.331*0 


where  SE  = 0.527  and  r = O.I+982 


3.  Step  J. — All  variables  were  then  moved  to  the  right  side,  and  a 
third  regression  was  performed  to  adjust  the  exponents  and  the  slope 
intercept : 


RD  = 0.366 


0.51*39  to.  531*3  r0.161+7 

P , 

.0.6911*  _0.  331*0 

t c2 


l.ll+io 


where  SE  = 0.520  and  r = 0.7820 


In  this  particular  version,  the  variable  P is  expressed  in  kips  (P„). 

K. 

1*.  Step  2a The  term  "log  t"  was  substituted  for  the  term  "t", 

and  Step  2 was  repeated: 


.0.5777 


RD 

0.51*39 


AW  ■ 00169 


(log  t) 


1.6676  ,0.3569 


where  SE  = O.526  and  r = 0.5035 
5.  Step  3a.  Step  3 was  repeated: 


RD  = 0.00062 


1 p \o. 51*39 


t0.5777  R0. 161*7 


(log  t)1'6676  cO-3569 


1.11*56 


where  SE  = 0.518  and  r = 0.781*2 


Note:  A higher  r value  was  attained  for  Steps  2a  and  3a  than  for 
Steps  2 and  3;  therefore,  the  imposed  term  "log  t"  was  allowed  to  remain. 
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6.  Step  k. — The  variable  P in  pounds  was  changed  to  in 

kips  to  make  the  coefficient  (0.00062)  higher: 


RD  = 0.01+59 


[Pk\0.5U39  t°'5777  E°-16hl 

fc)  (lo/t )l-«76  C?'35^ 


1.1*»57 


where  SE  = 0.51T  and  r = 0.781+2 


7.  Step  5. — The  exponents  of  variables  were  multiplied  by  1.11+57: 


RD  = 0.01+59 


, .0.6231 

K 

, 0. 6619  do.i887 

t r\ 

. p 

pi 

(log  t)1-*06  c"-1*089 

where  SE  = 0.517  and  r = 0.781+2 


Using  the  residuals  of  the  equation  from  Step  5,  1+5  points  that 

were  considered  undesirable  were  deleted  resulting  in  25l+  data  points. 

The  criteria  in  choosing  the  equation  from  Step  5 as  the  best  of  the 

equations  derived  were  first  that  this  equation  showed  an  increase  in 

RD  with  an  increase  in  P , t , and  R and  a decrease  in  RD  with 

P 

an  increase  in  , and  t as  should  be  expected.  In  addition, 

the  exponents  of  the  variables  seemed  reasonable  with  the  possible 
exception  of  the  exponent  for  repetitions  R , which  seemed  to  be  low. 
However,  this  equation  was  considered  the  best  overall. 

After  the  1+5  points  were  deleted,  a new  set  of  equations  was  de- 
rived from  the  remaining  25I+  data  points.  The  equation  considered  the 
best  of  these  was 


1+0 


0. 281+8 

R0.2UT6 

0020  0. 9335 

L1 

where  SE  = O.29I+  and  r = 0.9177 
or  by  moving  Cg  to  the  denominator 


0.1*707 


RD  = 0.171*1 


(log  t) 


2. 


RD  = 0.171+1 


p0. 1*707  t0.  5695  R0. 21*76 
(log  t)2-0020  c0.9335  <,0. 231*8 


Other  equations  that  were  considered  and  eliminated  were  as  follows 


RD  = 0.0723 


p0.  51*01  t0.  6666  R0. 2533 

C°’9303  (logPt)2-0U6^  C°-1U8 


where  SE  = 0.292  and  r = 0. 9187 
exponent  of  C 2 too  low 


RD 


0.121* 


p0.3933  .0.1*9 61*  0. 1829 
K \ R 

r O.8898  .0.6389  ^0.2728 
1 z c2 


where  SE  = 0.310  and  r = 0.9087 

derived  by  steps;  exponent  of  P and  R too  low 


1*1 


RD  = 1.075 


0.2328  0.3088  0. 2221 

rK  Lp  ^ 

^r5S59777^o33or 

where  SE  = 0.325  and  r = 0.8984 
exponent  of  P too  low 


RD  = 0.0348 


3o. 3987 

K 


0.4190 


"2/ 


C?'9293  (log  t) 


r0.2377 

1.8823 


where  SE  = 0.302  and  r = O.9125 

exponents  of  P and  t /C_  too  low 

P 2 


RD  = 0.156  -£ 

Cl, 


0.6942,t  ,0.3724 

k 


R 


0.2082 


(log  t)1'8888 


where  SE  = 0.334  and  r = 0.8921 


exponent  of  R too  low,  r lower  than  others 


RD  = O.OO676 


30.4706 

K 


0.2848 


,0.2476 


c0. 9335  (log  t)2.0019 


where  SE  = 0.294  and  r = 0.9176 

by  changing  P to  P^  final  equation  derived 


,0.1+483 


RD  = 1.510 


0.26l6 


,0.2548 


-0.9450  ,0.8125 
L1  1 


where  SE  = 0.297  and  r = 0.9159 


poor  results  with  substitution  of  "t"  instead  of 


'log  t' 


RD  = 0.0186 


po. 5170  to. 6191  R0. 2607 

E 

-0.9502  0.8311  -0.1216 
L1  t C2 


where  SE  = 0.295  and  r = 0.9170 


exponent  of  Cp  too  low. 


Two-Layer  Flexible  Pavement  Model 


The  data  search  resulted  in  the  selection  of  630  data  points  for 
flexible  pavements.  The  data  are  tabulated  in  Appendix  III.  Two 
hundred  and  ninety-three  of  the  data  points  were  applicable  to  three- 
layer  flexible  pavements. 

The  remaining  337  points  represented  a two-layer  flexible  pave- 
ment, or  one  not  having  a subbase  directly  above  the  subgrade.  The 
variables  considered  in  development  of  the  two-layer  flexible  pavement 
model  are  shown  in  Table  2.  Regression  analysis  was  performed  on  the 
337  data  points  in  steps  as  previously  described.  Several  forms  were 


imposed  upon  some  of  the  variables,  bub  certain  of  the  data  appeared  to 
be  unreasonable  in  terms  of  conventional  pavement  design.  The  residual 
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TABLE  2.— VARIABLES  FOR  TWO-LAYER  FLEXIBLE 
PAVEMENT  MODEL 


Index 

1 

2 

3 

4 

5 

6 

7 

8 


Variable 

RD  = rut  depth,  in. 

P = equivalent  single-wheel  load 
( ESWL ) , lb 


t 

t 

t 

C 

c 


= tire  pressure,  psi 
^ = AC  thickness,  in. 

2 = thickness  of  base,  in. 

^ = CBR  on  top  of  base 

2 = CBR  on  top  of  subgrade 

R = repetitions  of  load  or  passes 


error  values  between  these  preliminary  equations  and  the  actual  data 
were  used  to  eliminate  some  of  these  points.  A total  of  60  points  were 
eliminated  in  this  manner,  and  after  one  additional  iteration,  one  more 
data  point  that  had  been  overlooked  was  removed.  The  final  model  was 
therefore  based  upon  27 6 data  points.  Some  of  the  equations  developed 
and  rejected,  along  with  reasons  for  rejecting  them  are  shown  as  follows: 


RD  = 0.000831 


P 


pl.fcOO  r0.319 

0.306  ,0.577  r 1.206 
tl  2 °1 


0.17^ 

2 


where  SE  = 0.1*03  and  r = 0.8834 

t in  denominator,  indicating  that  rut  depth 

decreases  as  t increases,  which  is  erroneous. 
P 


1+4 


RD  = 0.00250 


.0.230 

tl 


P1 -269 

.0.755 

l2 


R°-319 

1.756  p0.0^0 
4 L2 


where  SE  = 0.1+21  and  r = 0.8720 


omitted  t ; exponent  of  was  low. 

... 


RD  = 0.003^1 


4.285 


+ t2) 


1.027 


R°-321 

1.700 

4 


0.077 

2 


where  SE  = 0.1+16  and  r = 0.8752 

same  as  above 


RD  = 0.001*30 

(l.25t1 


pl. 276  r0.320 
+ ^jl.033  ,1.703 


0. 072 
L2 


where  SE  = O.U15  and  r = 0.8755 


higher  r than  with  (t^  + t2)  but  exponent  of 
still  low 


These  examples,  along  with  several  others,  represent  the  efforts  made  to 
develop  the  correlation.  The  significant  factor  in  this  development, 
however,  was  combining  the  thicknesses  of  the  surface  and  base  course 
layers.  The  final  selected  model  uses  the  combined  form  of  the  thickness 
and  indicates  a common  log  relationship.  Additionally,  the  sensitivity 
of  the  subgrade  CBR  was  increased  to  a more  acceptable  level. 


*»5 


The  final  form  of  the  two-layer  flexible  pavement  model  is  given 


as  follows : 


RD  = 1.9431 


pl. 3127  t0.0499  r0.324o 

h e_ 


[log  (i.25tx  + t2; 


,3.4204  1.6877 


,0.1156 

2 


where  SE  = O.Ull  and  r = O.8779. 

This  version  represents  the  best  combination  of  variable  co- 
efficients and  apparent  behavior. 


Three-Layer  Flexible  Pavement  Model 

The  variables  pertinent  to  the  three-layer  flexible  pavement 
model  are  shown  in  Table  3.  The  three-layer  model  includes  nine 
independent  variables,  whereas  the  two-layer  model  only  included  seven 

TABLE  3.— VARIABLES  FOR  THREE-LAYER  FLEXIBLE  PAVEMENT  MODEL 


Index 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Variable 

RD  = rut  depth,  in. 

P = equivalent  single- wheel  load 
(ESWL),  lb 

t = tire  pressure,  psi 
t = thickness  of  AC,  in. 
t0  = thickness  of  base,  in. 

= CBR  on  top  of  base 
t^  = thickness  of  subbase,  in. 

C0  = CBR  on  top  of  subbase 
C^  = CBR  on  top  of  subgrade 
R = repetitions  of  load  or  passes 
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independent  variables.  The  two  additional  variables  are  subbase  thick- 


ness (t^)  and  subbase  CBR  (Cg),  respectively. 

Included  in  the  three-layer  flexible  pavement  data  was  one  test 
item,  which  was  a multiple-layered  system  (more  than  three  layers).  For 
the  data  analysis,  the  values  used  for  the  base,  subbase,  and  subgrade 
CBR  variables  of  the  62  data  points  from  the  multiple-layered  item  were 
usually  an  average  CBR  of  two  layers.  Similarly,  the  thickness  variables 
were  usually  the  sum  of  the  thicknesses  of  the  two  layers.  Since  these 
values  may  not  have  accurately  represented  the  actual  pavement  structure 
of  the  test  item,  the  62  data  points  were  removed  to  see  if  the  corre- 
lation coefficient  would  improve.  A better  correlation  was  obtained,  and 
these  62  points  were  eliminated  from  further  analysis. 

As  in  the  case  of  the  previous  equations,  numerous  computer  runs 
were  tried  using  the  remaining  231  points  to  obtain  an  equation  con- 
sidered the  best  fit.  Some  of  the  equations  generated  and  rejected, 
along  with  reasons  for  rejection,  are  as  follows: 


RD  = 0.000021 


pl. 158  t0.256  R0. 376  t0.503 

.0.53^  p0.  39^  . 0.  S27  0.5M  p0. 289 
2 1 2 2 3 


where  SE  = 0.528  and  r = 0.8217 


t^  in  the  numerator  indicating  improper  relation 
to  RD 


i<7 


RD  = 0. 000001 T 


.,1.625  ,0.1l4  0. 367 


0.531  ,,0.682  0.590  .0.595  n0.203 

2 L1  S C2  C3 


3 ~2 

where  SE  = 0.552  and  r = 0.8037 
omitted  t^  in  attempt  to  improve  correlation 


RD  = 0.000011 


1.506  /log  t \0*291 


1 + log  t 


R 


0.330 


II  ' 1_  ' 

t0.4l0  t0.523  C0.617  0.164 


where  SE  = 0.635  and  r = 0.7281 
exponent  of  too  low,  r lower. 


RD  = 0.000005 


p I1-516 //si 

0.276 

D0.  330 

l1  + log hi  Is  i 

R 

where  SE  = O.638  and  r = 0.7258 
r still  low 


The  final  equation  selected  was  one  of  three  equations  considered  almost 
equally  good.  To  obtain  a higher  coefficient,  the  ESWL  (P)  was  changed 
from  pounds  to  kips  (P^)  in  each  of  the  following  three  equations: 
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RD  = 0.2218 


RD  = 0.2181* 

where  SE  = O.566  and  r = 0.7920 


RD  = 0.2065 

where  SE  = O.568  and  r = 0.7907 


In  the  first  of  the  three  equations  above,  it  can  be  noted  that  for  two 
different  pavement  structures,  as  long  as  the  sum  of  t^  and  t 2 are 
equal  it  does  not  matter  what  the  individual  thicknesses  of  t^  and  t^ 
are;  thus,  no  benefit  is  to  be  gained  by  using  asphaltic  concrete  in  lieu 
of  base  course  material.  This  was  considered  erroneous.  However,  in  the 
other  two  equations,  if  the  sum  of  t^  and  t 2 are  equal,  the  pavement 
with  the  larger  asphalt  thickness  would  decrease  the  rut  depth.  Since 
the  second  equation  has  a higher  correlation  than  the  third,  it  was 
selected  as  the  best  equation.  Using  the  residuals  of  the  second  equa- 
tion, 25  points  that  were  considered  erroneous  were  deleted,  resulting  in 
206  data  points. 

1*9 


„i. 61U  .0.091+  D0. 363 

PK  tp  R 

(1.25tl  * t,)0-257  cj-577  t»'523  C2'545 


,0.337 


^1 . 638  0.085  D0. 363 

r l r 

K p 


<1.25^  ♦ t2)°'335  C?'8* 


,0. 302 


f 


The  selection  of  the  best  equation  using  the  remaining  206  points 
was  again  narrowed  to  one  of  three  equations.  The  first  equation  was  the 
same  form  as  that  selected  earlier. 


RD  = O.I8U6 


(l.25t1 


+ t2) 


1.709 

K 

o^of~c 


t0.050 

_E_ 

0.761 
1 t 


ro.3i+6 

0.582 

3 


0.571 

2 


O.269 

3 


where  SE  = 0. 1+1+5  and  r = 0.81+lU 


Substituting  the  common  log  of  the  thicknesses  for  the  thickness  vari- 
ables resulted  in  the  other  two  equations  cons.i’^red  best. 


RD 


o.o68op^,T3° 
[log  (l.25t1  + t2)]1,223 


t0.0l+5  R0.  3^6 

, 

,-,0. 759  . 0.  589  „0.  578 


O.2J+3 

3 


where  SE  = 0. 1+1+3  and  r = 0.81+25 


0.0312P; 


RD  = 


1.525 

K 


t0.090  R0.3l+5 

2 

1.167  c0.551  c0.309 


°-465  doe  t 


[log  (l.25t  + 


t2)] 


where  SE  = 0. 1+1+1+  and  r = 0.81+18 

All  three  of  the  above-mentioned  equations  showed  an  increase  in  rut 
depth  proportional  to  an  increase  in  P , t , and  R and  a decrease 
in  rut  depth  proportional  to  an  increase  in  the  other  six  variables.  The 
final  selection  of  the  best  of  these  three  models  was  delayed  until  the 
analysis  of  the  other  group  of  flexible  pavement  data  was  performed.  The 
latter  of  these  three  models  was  then  selected  based  upon  the  behavior  of 
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the  thickness  variables  as  well  as  its  similarity  to  the  two-layer 
flexible  pavement  model.  Carrying  the  exponents  to  four  decimal  places 
and  rearranging,  the  final  equation  becomes 


RD  - 0.03117  |pJ'5255  t°'0®97  R9-3“50/Uo8  (1.25t1+  t2)]0'881*7 

* (log  C°-7616  a"'5505  C°'3089} 

where  SE  = 0 . UUU  and  r = O.8U18 


DEVELOPMENT  OF  RELIABILITY  MODELS 


Development  Procedure 

The  deterioration  models  for  rut  depth  analyses  of  the  four  respec- 
tive pavement  types  as  developed  previously  are  necessarily  deterministic 
models.  They  represent  not  necessarily  the  best  fit  of  all  the  data,  but 
instead  a "good  fit"  as  indicated  by  the  error  and  correlation  values. 
Therefore,  the  models  predict  rutting  in  a sense  that  the  predictive 
error  could  be  either  positive  or  negative  depending  upon  the  specific 
set  of  data  being  used.  To  further  expand  the  applicability  and  utility 
of  the  models,  statistical  reliability  concepts  are  invoked  to  account 
for  the  variability  of  all  input  data  and  predict  rutting  in  terms  of 
expected  rut  depth  and  variance  from  that  expected  rut  depth.  This  con- 
cept provides  for  models  for  each  pavement  type  that  utilizes  a total 
description  of  the  input  variables  in  terms  of  means  and  variances  of  a 
set  of  values  for  one  variable.  This  further  accounts  for  the  variabil- 
ity of  pavement  properties  in  a statistical  manner  instead  of  accept- 
ing only  averages  of  values  of  measured  parameters,  such  as  thickness 
and  CBR.  The  statistical  determination  of  the  rut  depth  in  terms  of 
expected  value  and  variance  provides  not  only  for  better  analysis  of  rut 
depth  in  terms  of  an  expected  rut  depth  and  a probable  deviation  from 
that  value  but  also  for  an  accounting  for  material  variability  in  a 
statistical  manner.  The  primary  benefit  of  using  total  input  data  and 
determining  the  expected  values  and  variances  of  a probability  density 


function  of  rut  depth  is  the  ability  to  evaluate  the  reliability  of 
the  pavement.  The  reliability  of  a facility  can  be  determined  in 
terms  of  the  probability  that  a given  rut  depth  will  occur  under  given 
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circumstances . It  follows  then  that  such  a system  can  he  an  excellent 
evaluation  tool  as  well  as  a design  tool.  The  great  advantage  as  a de- 
sign procedure  is  the  capability  to  adjust  reliability  or  conservatism 
to  a desired  value  and  to  select  design  parameters  to  suit  those  con- 
ditions. To  address  these  concepts  of  reliability  and  the  accounting 
for  material  variability,  it  was  necessary  to  develop  stochastic  models 
for  the  definition  of  a probability  density  function  of  the  rut  depth. 

Probabilistic  pavement  design  concepts  have  been  applied  to  pave- 
ment studies  since  the  1960's.  References  5,  7,  18,  21,  23,  and  2k  give 
illustrations  of  these  studies  wherein  the  concept  of  pavement  relia- 
bility has  been  adopted.  The  reliability  of  a pavement  is  a statis- 
tical measure  of  the  probability  that  a pavement  will  perform  in  a given 
manner  during  its  life.  Lu,  Lytton,  and  Moore  (20)  showed  the  use  of  the 
above-mentioned  expected  value  and  variance  equations  as  applied  to  fore- 
casting serviceability  loss  in  pavements.  These  principles  were  adopted 
and  utilized  to  develop  expected  value  and  variance  equations  for  pre- 
dicting change  in  rut  depth.  The  methods  for  the  determination  of 
expected  value  and  variance  of  rut  depth  by  a Taylor's  series  expansion 
(20)  are  shown  in  the  following  paragraphs. 

If  rut  depth  (RD)  is  considered  a continuous,  random  variable  with 
some  probability  density  function  f(RD)  , the  expected  value  of  H(RD) 
is  defined  as 

00 

E[H(RD) ] = f H(RD)  f (RD)  dRD 
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The  expected  value  of  RD  is  the  mean  or  average  of  RD  , which  is 

termed  uRR  , or  it  can  be  written  as  uRD  = E(RD)  . The  variance  of  RD 

2 

as  a variable  is  denoted  by  and  is  defined  to  be 

KD 

°RD  = E^RD  - J 

Its  positive  square  root  is  the  standard  deviation  of  RD  . Thus, 

°RD  = i°KD 

The  operation  of  taking  expected  values  and  variances  of  random  variables 
is  found  in  various  textbooks  and  is  illustrated  by  Lu,  Lytton,  and 
Moore  (20).  Occasionally,  taking  the  expected  value  of  a complicated 
function  can  be  a difficult  process  as  in  the  case  of  the  rut  depth 
models  previously  shown.  In  order  to  overcome  these  difficulties,  the 
expected  value  was  approximated  by  taking  a Taylor's  series  expansion  and 
truncating  all  but  the  first  three  terms  as  follows: 

f(RD)  = f(RD  - ARD)  + f ' (RD  - ARD)ARD  + 1/2  F"(RD  - ARD)ARD2  + ... 

If  ARD  becomes  RD  - pRR  , then  the  final  generalized  form  can  be  ex- 
pressed as  follows: 

E[f(RD)  I . f(„BD)  ♦ 1/2  f"(uRD)<4 

Taking  the  variance  of  the  rut  depth  models  was  also  a painstaking 
operation,  further  complicated  by  the  forms  imposed  upon  some  of  the 


variables.  The  Taylor's  series  expansion  was  again  applied.  As  was 


previously  stated,  the  variance  was 


°RD  " ^RD  ^RD^  ^ 


The  variance  of  the  rut  depth  models  is  denoted  V(RD)  and  is  expressed 


as 


V(RD)  [f'^RD^  °RD  “ T7  f1  (pRc^  °RD  + °J!.of 


where  o.  . is  the  variance  of  lack  of  fit. 

Hof 

The  deterministic  rut  depth  analysis  models  have  already  been  de- 
rived in  this  document.  The  principles  shown  above  were  invoked  to  de- 
velop the  expected  value  and  variance  models  from  the  original  model. 

The  nature  of  the  equations  that  provided  for  numerous  mathematical  oper- 
ations precludes  showing  derails  of  development  in  this  document.  How- 
ever, these  mathematical  operations  are  considered  fundamental  and  are 
not  shown.  The  final  forms  of  the  equations  are  listed  in  the  following 
figures. 


Expected  Value  and  Variance  Models 

Figs,  j.6-19  show  the  expected  value  and  variance  models  for  the 
four  respective  categories.  The  rut  depth  models  previously  discussed 
are  also  shown.  This  version  of  the  rut  depth  model  differs  from  the 
original  model  in  that  each  independent  variable  actually  represents 
a mean  value,  and  the  equation  is  therefore  redesignated  "Q." 
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RD  = 0.1101 


pO.4925  t^8340  (log  t)0A293  R°-5018 


-1.9773  -1.2015 
C2 


(=  Q in  the  following  equations) 
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(=  Q in  the  following  equations) 
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FIG.  IT-— GRAVEL-SURFACED  FACILITY  EXPECTED  VALUE, 
VARIANCE,  AND  RUTTING  MODELS 
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RD  = 1.9431 
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FIG.  18.— TWO-LAYER  FLEXIBLE  PAVEMENT  EXPECTED  VALUE, 
VARIANCE,  AND  RUTTING  MODELS 
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RD  = 0.03117  jpl-5255t0.°897Rc|.31450^ + t^]0'88**7 

« (log  tj)1'167*1  cj-7616  C°'5505  C°'30(4 


(=  Q in  the  following  equations) 


FIG.  19.— THREE- LAYER  FLEXIBLE  PAVEMENT  EXPECTED  VALUE, 
VARIANCE,  AND  RUTTING  MODELS 
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C\J  4-> 


Actually  "Q"  exists  as  a factored  term  from  the  expected  value  and 
variance  models,  developed  in  the  process  of  evaluating  the  first  and 
second  partial  derivatives  of  the  original  model. 

As  can  he  seen,  the  models  are  combined  equations  representing 
the  first  three  terms  of  a Taylor's  series  expansion.  Many  of  the 
terms  have  been  allowed  to  remain  in  unfactored  form  to  display  the 
actual  nature  of  the  models  to  the  reader.  Any  attempt  to  solve  the 
models  by  hand  should  be  preceded  by  as  many  simplifications  of  the 
models  as  possible.  However,  hand  solution  of  the  models  is  necessarily 
time-consuming  and  leaves  room  for  many  possible  errors.  The  models  have 
been  programmed  for  computer  solution  as  part  of  the  analysis  model,  as 
shown  in  Appendix  IV. 

Determination  of  Reliability 

It  should  be  noted  that  each  variable  in  the  two  models  is  used 
in  terms  of  its  respective  mean  and  variance  and  that  in  turn  an  expected 
value  and  variance  of  rut  depth  is  determined.  As  has  been  previously 
stated,  reliability  as  defined  and  used  in  this  study  is  the  probability 
that  the  rut  depth  will  not  exceed  some  predetermined  value  subject  to 
conditions  that  are  expressed  by  the  independent  variables. 

The  nature  of  the  data  including  the  range  and  distribution  will 
be  discussed  in  succeeding  parts  of  this  document,  along  with  an  explana- 
tion for  assumption  of  normal  distribution  on  the  dependent  variable.  If 
normality  is  assumed,  the  reliability  statistic  P used  to  determine  the 
reliability  R that  the  rut  depth  will  not  exceed  some  maximum  value 
RD^  can  be  expressed  by 
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where 


RD  - E(RD) 

P = ' 

/v(rd) 


E(RD)  = the  expected  value  of  the  rut  depth  as  determined  from  the 
appropriate  model 

V(RD)  = the  variance  of  the  rut  depth  as  similarly  determined 
Figure  20  illustrates  a standard  normal  distribution  curve  and  the  param- 
eters utilized  to  compute  the  reliability  statistic  P . Table  b is  a 
table  of  areas  under  the  standard  normal  distribution,  or  its  entries  are 
the  values  of  reliability  R . 


PVVfRD) 

E(  RD)  +E(RD) 

h-* -t- 

FIG.  20.  STANDARD  NORMAL  DISTRIBUTION  CURVE  ILLUSTRATING 
E(RD)  , /vTrdT  , and  RD^ 

In  order  to  determine  reliability  R from  Table  4 or,  that  is,  the  area 
under  the  distribution  curve  defined  by  E(RD)  and  V(RD)  and  to  the 
left  of  maximum  rut  depth  RDA  , enter  the  table  with  the  value  of  P 
determined  previously  and  read  the  area  under  the  distribution  curve  that 
is  the  reliability.  As  an  illustrative  example,  assume  the  following 
values: 
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TABLE  It.— NORMAL  DISTRIBUTION  FUNCTION 
F(P)  - R 


p 

0.00 

0.01 

0.02 

0.03 

0.0U 

0.05 

0.06 

0.07 

0.08 

0.09 

0.0 

0.5000 

0.501*0 

0.5080 

0.5120 

0.5160 

0.5199 

0.5239 

0.5279 

0.5319 

0.5359 

0.1 

0.5398 

0.5*t38 

0.51*78 

0.5517 

0.5557 

0.5596 

0.5636 

0.5675 

0.5711* 

0.5753 

0.2 

0.5793 

0.5832 

0.5871 

0.5910 

0.591*8 

0.5987 

0.6026 

0.6061* 

0.6103 

0.6ll*l 

0.3 

0.6179 

0.6217 

0.6255 

0.6293 

0.6331 

0.6368 

0.61*06 

0.61*1*3 

0.61*80 

0.6517 

O.U 

0.655** 

0.6591 

0.6628 

0 . 666U 

0.6700 

0.6736 

0.6T72 

0.6808 

0.68U1* 

0.6879 

0.5 

0.6915 

0.6950 

0.6985 

0.7019 

0.7051* 

0.7088 

0.7123 

0.7157 

0.7190 

0.7221* 

0.6 

0.7257 

0.7291 

0.732** 

0.7357 

0.7389 

0.71*22 

0.71*51* 

0.71*86 

0.7517 

0.751*9 

0.7 

0.7580 

0.7611 

0.761*2 

0.7673 

0.7701* 

0.773>* 

0.7761* 

0.7791* 

0.7823 

0.7852 

0.8 

0.7881 

0.7910 

0.7939 

0.7967 

0.7995 

0.8023 

0.8051 

0.8078 

0.8106 

0.8133 

0.9 

0.8159 

0.8186 

0.8212 

0.8238 

0.8261* 

0.8289 

0.8315 

0.831*0 

0.8365 

0.8389 

1.0 

0.81*13 

0.81*38 

0.81*61 

0.81*85 

0.8508 

0.8531 

0.8551* 

0.8577 

0.8599 

0.8621 

1.1 

0.861*3 

0.8665 

0.8686 

0.8708 

0 . 8729 

0.871*9 

0.8770 

0.8790 

0.8810 

0.8830 

1.2 

0.881*9 

0.8869 

0.8888 

0.8907 

0.8925 

0.891*1* 

0.8962 

0.8980 

0.8997 

0.9015 

1.3 

0.9032 

0.901*9 

0.9066 

0.9082 

0.9099 

0-9U5 

0.9131 

0.911*7 

0.9162 

0.917T 

1.1* 

0.9192 

0.9207 

0.9222 

0.9236 

0.9251 

0.9265 

0.9279 

0.9292 

0.9306 

0.9319 

1.5 

0.9332 

0.931*5 

0.9357 

0.9370 

0.9382 

0.939>* 

0.91*06 

0.91*18 

0.91*29 

0.91*1*1 

1.6 

0.91*52 

0.91*63 

0.9l*7l* 

0.91*81* 

0.9U95 

0.9505 

0.9515 

0.9525 

0.9535 

0.951*5 

1.7 

0.955** 

0.956** 

0.9573 

0.9582 

0.9591 

0.9599 

0.9608 

0.9616 

0.9625 

0.9633 

1.8 

0.961*1 

0.961*9 

0.9656 

O.966U 

0.9671 

0.9678 

0.9686 

0.9693 

0.9699 

0.9706 

1.9 

0.9713 

0.9719 

0.9726 

0.9732 

0.9738 

0.971*1* 

0.9750 

0.9756 

0.9761 

0.9767 

2.0 

0.9772 

0.9778 

0.9783 

0.9788 

0.9793 

0.9798 

0.9803 

0.9808 

0.9812 

0.9817 

2.1 

0.9821 

0.9826 

0.9830 

O.983I* 

O.9838 

0.981*2 

0.981*6 

0.9850 

O.985I* 

0.9857 

2.2 

0.9861 

O.986U 

0.9868 

0.9871 

0.9875 

0.9878 

0.9881 

O.988I* 

0.9887 

0.9890 

2.3 

0.9893 

0.9896 

0.9898 

0.9901 

0.9901* 

0.9906 

0.9909 

0.9911 

0.9913 

0.9916 

2.1* 

0.9918 

0.9920 

0.9922 

0.9925 

0.9927 

0.9929 

0.9931 

0.9932 

0.9931* 

0.9936 

2.5 

0.9938 

0.991*0 

0.991*1 

0.991*3 

0.991*5 

0.991*6 

0.991*8 

0.991*9 

0.9951 

0.9952 

2.6 

0.9953 

0.9955 

0.9956 

0.9957 

0.9959 

0.9960 

0.9961 

0.9962 

0.9963 

0.9961* 

2.7 

0.9965 

0.9966 

0.9967 

0.9968 

0.9969 

0.9970 

0.9971 

0.9972 

0.9973 

0.9971* 

2.8 

0.9971* 

0.9975 

0.9976 

0.9977 

0.9977 

0.9978 

0.9979 

0.9979 

0.9980 

0.9981 

2.9 

0.9981 

0.9982 

0.9982 

0.9983 

O.998I* 

0.9981* 

0.9985 

0.9985 

0.9986 

0.9986 

3.0 

0.9987 

0.9987 

0.9987 

0.9988 

0.9988 

0.9989 

0.9989 

0.9989 

0.9990 

0.9990 

3.0 

0.9990 

0.9991 

0.9991 

0.9991 

0.9992 

0.9992 

0.9992 

0.9992 

0.9993 

0.9993 

3.2 

0.9993 

0.9993 

0.999** 

0 . 9991* 

0.9991* 

0.9991* 

0.999!* 

0.9995 

0.9995 

0.9995 

3.3 

0.9995 

0.9995 

0.9995 

0.9996 

0.9996 

0.9996 

0.9996 

0.9996 

0.9996 

0.9997 

3.U 

0.9997 

0.9997 

0.9997 

0.9997 

0.9997 

0.9997 

0.9997 

0 . 9957 

0.9997 

0.9998 

6U 
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RDa  = 3 in. 

E(RD)  = 2 in. 

V(RD)  = 1 in. 

RDa  - E(RD) 

then  P = = = — = 1.  If  Table  4 is  entered  with  a 

/V(RD)  JT 

value  of  1 in  the  left-hand  column,  interpolation  is  not  necessary. 
Also,  it  can  immediately  be  seen  that  R = 0.8413,  which  means  that 
there  is  a probability  of  0.8413  that  the  rut  depth  will  not  exceed  a 
predetermined  maximum  value  of  3 in.  Table  4 is  typical  of  similar 
tables  found  in  most  statistics  textbooks.  Details  for  the  operation 
of  the  computer  programs  are  presented  in  Appendix  IV. 
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DIFFERENTIAL  ANALYSIS  SYSTEM 

Description 

Fig.  21  gives  a logic  diagram  of  DAS  as  defined  in  this  document. 

The  system  provides  for  a utilization  of  the  rutting  models  to  determine 
the  rate  of  deterioration  and/or  reliability  of  any  of  the  four  types  of 
facilities  described  in  terms  of  rutting.  The  term  "differential"  has 
been  given  to  the  computational  system  to  emphasize  the  fact  that  dif- 
ferences in  results  caused  by  changes  in  input  can  be  determined  by  the 
user  in  any  assessment  of  damage  caused  by  various  vehicle  types  or  the 
effects  of  changes  in  the  structure.  The  DAS,  as  shown  in  Fig.  21,  pro- 
vides for  one  automatic  iteration  of  the  computational  processes.  Dif- 
ferential analysis  as  described  above  can  be  achieved  simply  by  repeated 
iterations  of  the  system  while  changing  any  variable  or  variables  de- 
sired. The  system  as  shown  is  adequate  for  limited  use  where  the  various 
models  apply  and  is  adequate  to  develop  the  original  hypothesis  that 
life-cycle  management  can  be  achieved  througn  deterioration  and  reli- 
ability concepts. 

The  DAS  is  programmed  for  computer  solution.  The  program  listing, 
input  listing,  output  listing,  instructions  for  use,  and  some  example 
problems  are  given  in  Appendix  IV.  Appendix  IV  should  be  utilized  by  any 
reader  concerned  with  operating  the  system.  The  following  is  a descrip- 
tion of  the  DAS  (Fig.  21).  Block  numbers  identify  the  particular  block 
in  the  logic  diagram  being  discussed. 

1.  Blocks  1 and  2. — Blocks  1 and  2 provide  for  input  of  data 
describing  each  independent  variable  considered.  The  data  can  be 
entered  in  deterministic  single-value  form  but  should  be  entered  in  terms 
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FIG.  21. —DIFFERENTIAL  ANALYSIS  SYSTEM  (DAS),  LOGIC  DIAGRAM  (STAGE  l) 
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of  the  mean  and  variance  of  a group  of  data  in  order  to  fully  utilize  the 
potential  of  the  DAS. 

2.  Block  3. — Block  3 consists  of  the  models  developed  earlier  in 
this  report,  including  the  rutting,  variance,  and  expected  value  models. 
The  models  for  all  four  types  of  structures  are  included,  and  the  selec- 
tion of  the  appropriate  set  is  incumbent  to  using  the  DAS.  Block  3 is 
the  point  where  RD  , V(RD)  , and  E(RD)  , respectively,  are  computed. 

3.  Block  4. — Block  ^ provides  for  input  of  the  allowable  rut  depth 
and  an  allowable  (minimum)  reliability  level.  The  actual  input  is  at  the 
beginning  of  the  system  but  is  shown  here  as  its  first  point  of  use. 

U.  Block  5 . — Data  are  available  at  this  point  to  determine  the 
reliability  of  the  facility.  This  is  accomplished  as  shown  in  Block  5 
and  in  the  earlier  paragraph  on  reliability  determination. 

5.  Block  6. — If  the  problem  is  a differential  analysis  problem, 
Block  6 is  the  point  where  the  decision  is  made  to  perform  additional 
iterations  and  obtain  new  rut  depth  and  reliability  values  that  reflect 
effects  of  changes  made  in  any  variables. 

6.  Blocks  7 and  8. — These  blocks  provide  for  change  in  input  to  the 
DAS  in  accordance  with  the  type  of  differential  analysis  being  conducted. 

7.  Block  9. — Block  9 is  an  additional  decision  point  that  provides 
for  new  iterations  when  the  limiting  conditions  of  either  rut  depth  or 
reliability  have  not  been  met. 

8.  Block  10. — Block  10  provides  for  an  output  listing  of  the  data 
obtained.  These  data  are  actually  accessed  at  the  decision  points  to 
determine  need  and  desirability  of  performing  additional  iterations. 

Fig.  22  shows  the  logic  of  an  expanded  differential  analysis  system 

that  considers  surface  loss,  roughness,  and  cracking  as  well  as  rutting. 
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It  uses  the  same  logic  as  the  system  shown  in  Fig.  21  in  determining  de- 
terioration and  reliability  but  includes  traffic  design  and  economic 
considerations  also.  Consequently,  it  will  require  more  iterations  to 
determine  the  optimal  feasible  solution.  Although  not  considered  es- 
sential to  this  hypothesis  development,  it  is  shown  here  to  give  the 
reader  a concept  of  the  validated  and  automated  version  of  the  DAS  that 
will  only  be  utilized  by  other  agencies  as  the  DAS  is  adopted. 

Fig.  23  is  another  version  of  the  second  stage  of  the  DAS,  showing 


FIG.  23.— DIFFERENTIAL  ANALYSIS  SYSTEM.  SUMMARY 
DIAGRAM  OF  STAGE  2 


the  essential  elements.  The  system  will  provide  for  optimization  of 
-itt ing  and  reliability  values  with  respect  to  constraints  imposed 

.►'h  the  input  of  values  or  ranges  of  values  representing  the  vari- 
,-iiitional  optimization  will  be  achievable  in  terms  of  design 
• it ion  costs  as  well  as  selection  among  the  four  pavements 
A a if. a pertinent  to  other  pavement  types,  especially  rigid 
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pavements  and  selected  hybrid  pavements,  become  available,  models  will 
be  included  for  their  analysis. 


Range  and  Distribution  of  Variables 

The  DAS  is  considered  to  be  applicable  to  design  and  evaluation 
problems  where  normally  encountered  values  of  the  variables  are  utilized. 
The  range  of  applicability  of  any  computational  system  which  uses  empir- 
ical formulas  as  developed  herein  is  constrained  by,  if  not  limited  to, 
the  range  of  observed  data  upon  which  it  is  based. 

Figs.  24-55  graphically  portray  the  range  and  distribution  of  vari- 
ables upon  which  the  DAS  is  based.  As  can  be  seen,  the  boundaries  will 
normally  encompass  most  vehicle  and  road  characteristics  to  be  input  as 
variables. 

The  data  are  not  normally  distributed  in  most  cases  as  was  indicated 
in  tests  for  normality.  However,  it  was  realized  that  the  data  collected 
and  utilized  were  not  intended  to  group  about  any  particular  value  as  a 
whole.  Instead,  the  data  shown  in  Appendix  I group  about  various  values 
within  each  given  series  of  tests  in  a normal  manner.  For  this  reason, 
normal  distribution  is  assumed  in  determining  reliability  as  has  pre- 
viously been  discussed.  Validation  of  the  DAS  for  use  by  various 
agencies  will  necessarily  include  tests  for  normality  or  determination 
of  the  nature  of  statistical  distribution  for  the  purpose  of  more  suc- 
cinctly determining  reliability.  In  general,  the  DAS  is  applicable 
within  data  ranges  and  especially  where  data  are  concentrated. 
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FIG.  27 . — UNSURFACED  FACILITY,  GRAVEL  THICKNESS  DISTRIBUTION 
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FIG.  29 . — UNSURFACED  FACILITY,  SUBGRADE  CBR  DISTRIBUTION 
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FIG.  32. --GRAVEL-SURFACED  FACILITY,  LOAD  DISTRIBUTION 


120  n 


100  H 


Base  Thickness  (t^),  in. 

FIG.  3^ . — GRAVEL-SURFACED  FACILITY,  GRAVEL  THICKNESS  DISTRIBUTION 
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FIG.  36.— GRAVEL-SURFACED  FACILITY,  SUBGRADE  CBR  DISTRIBUTION 


FIG.  37. --GRAVEL-SURFACED  FACILITY,  REPETITION  DISTRIBUTION 
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FIG.  38. — TWO-LAYER  FLEXIBLE  PAVEMENT,  RUT  DEPTH  DISTRIBUTION 


FIG.  1*1. — TWO- LAYER  FLEXIBLE  PAVEMENT,  ASPHALT  THICKNESS  DISTRIBUTION 
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FIG.  4 3. —TWO- LAYER  FLEXIBLE  PAVEMENT,  BASE  CBR  DISTRIBUTION 
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FIG.  bb . — TWO-LAYER  FLEXIBLE  PAVEMENT,  SUBGRALE  CBR  DISTRIBUTION 


FIG.  4 5.— TWO-LAYER  FLEXIBLE  PAVEMENT,  REPETITION  DISTRIBITION 


FIG.  k6. — THREE- LAYER  FLEXIBLE  PAVEMENT,  RUT  DEPTH  DISTRIBUTION 


FIG.  1*7. -THREE-LAYER  FLEXIBLE  PAVEMENT,  LOAD  DISTRIBUTION 
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FIG.  48. — THREE-LAYER  FLEXIBLE  PAVEMENT,  TIRE  PRESSURE  DISTRIBUTION 


FIG.  49. — THREE-LAYER  FLEXIBLE  PAVEMENT,  ASPHALT  THICKNESS  DISTRIBUTION 
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FIG.  52.— THREE- LAYER  FLEXIBLE  PAVEMENT,  SUBBASE  THICKNESS  DISTRIBUTION 


FIG.  53.— THREE- LAYER  FLEXIBLE  PAVEMENT,  SUBBASE  CBR  DISTRIBUTION 
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FIG.  5^ • — THREE-LAYER  FLEXIBLE  PAVEMENT,  SUBGRADE  CBR  DISTRIBUTION 


FIG.  55. --THREE- LAYER  FLEXIBLE  PAVEMENT,  REPETITION  DISTRIBUTION 
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Utilization  of  the  DAS 


There  are  numerous  applications  of  the  DAS  that  can  be  made  by  those 
concerned  with  life-cycle  management  or  any  aspect  thereof.  It  is  again 
pointed  out  that,  as  the  case  with  any  similar  system,  validation  is  re- 
quired to  render  the  DAS  directly  applicable  to  specific  locales  having 
unique  conditions.  Previously  discussed  data  collection  programs  that 
are  in  progress  are  expected  to  provide  data  for  this  validation  effort. 
The  data  collected  are  expected  to  provide  rate  of  change  in  rut  depth  as 
a function  of  vehicle  repetitions  under  several  conditions  of  road  types, 
traffic  levels,  and  climate.  These  rut  depths  will  be  predicted,  along 
with  reliability  levels,  using  the  DAS.  The  predicted  and  measured  rut 
depth  values  will  be  compared.  An  acceptable  comparison  will  constitute 
validation,  but  lack  of  acceptable  comparison  will  dictate  modification 
of  the  expected  value  and  variance  models  for  use  under  those  conditions. 
The  DAS  is,  however,  applicable  to  all  situations  in  a general  or  sto- 
chastic sense  and  can  be  used  accordingly  by  the  discerning  engineer. 

Since  this  version  is  in  terms  of  rut  depth  as  the  dependent  variable  and 
major  item  of  analysis,  then  rut  depth  is  the  primary  criterion  by  which 
to  judge  state  of  deterioration.  Validation  of  the  DAS  not  only  can 
achieve  local  applicability  but  can  incorporate  the  other  deterioration 
modes  as  a data  base  is  made  available.  Field  data  collection  for  vali- 
dation purposes  that  is  currently  under  way  has  been  discussed  previously. 

Although  the  specific  applications  of  the  DAS  are  numerous  and  nec- 
essarily depend  upon  user  needs,  some  of  the  more  prominent  applications 
are  described  as  follows: 

1.  Design  and  Evaluation.— The  DAS  is  applicable  to  CE  design  and 
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evaluation  problems  in  the  same  sense  as  are  current  criteria  due  to 
similarity  in  data  bases  and  results.  The  added  features  of  the  DAS  are 
namely  the  ability  to  (a)  modify  limiting  failure  criteria  (rut  depth), 
(b)  adjust  conservatism  to  any  desired  degree  by  imposing  a required 
degree  of  reliability,  and  (c)  determine  the  reliability  of  a facility. 

2.  Optimization. — Iterations  of  the  DAS  while  making  changes  in 
appropriate  variables  can  provide  for  optimization  of  a design  with  re- 
spect to  cost,  reliability,  serviceability,  layer  properties,  and 
materials. 

3.  Differential  Analysis. — Iterations  of  the  DAS  provide  directly 
for  the  analysis  of  the  effects  of  different  quantities  and  magnitudes 
of  loads.  The  equivalent  single-wheel  concept  (12)  makes  this  possible 
by  providing  the  capability  to  incorporate  various  vehicle  configura- 
tions. This  feature  provides  a quantitative  basis  for  assessment  of 
damages  caused  by  various  categories  of  vehicles  and,  when  used  on  a rel- 
ative basis,  would  not  require  locality  validation  of  the  DAS. 

U.  Planning. — The  DAS  can  be  considered  an  effective  stochastic- 
type  planning  tool  for  quantitative  estimation  of  future  maintenance  and 
repair  needs  as  well  as  time-to-maintenance  estimation.  This  feature,  in 
connection  with  such  procedures  as  Critical  Path  Method  (CPM)  and  Pro- 
gram Evaluation  and  Review  Technique  (PERT),  can  be  used  to  effectively 
program  work  loads  and  expenditures. 

5.  Military  Operations. — The  tactical  and  logistical  operations 
that  could  benefit  from  use  of  the  DAS  are  too  numerous  to  mention  in  de- 
tail. Such  considerations  as  optimization  of  construction  capabilities 
by  constructing  facilities  having  only  a required  reliability  and  using 
facility  reliability  concepts  to  aid  in  tactical  planning  and  maneuvers 
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are  key  considerations  that  could  be  better  quantified  using  DAS 

concepts.  ] 

The  use  of  the  DAS  as  it  exists  herein  necessarily  includes  use  of 
the  CBR  method  of  strength  evaluation,  which  is  in  itself  not  a true 
physical  material  parameter.  This  feature  is  not  to  be  considered  a 
deterrent,  however,  to  prospective  users  bound  to  other  design  pro- 
cedures. Material  strength  parameters  can  be  stated  in  any  suitable 
terms  where  a sufficient  data  base  exists  for  validation.  True  mate- 
rial parameters,  such  as  Poisson's  ratio  (p)  and  elastic  properties  j 

(E),  in  various  forms  can  be  utilized  and  would  provide  for  a more 
rational  approach  to  the  overall  operation.  Other  methods  of  portray- 
ing strength  can  be  used  where  data  are  available.  In  all  cases,  any 
bias  is  removed  in  the  actual  correlation  indicated  during  the  valida- 
tion stage  and  should  be  the  basis  of  judgment  as  to  whether  a partic- 
ular procedure  is  employed. 

Example  Problem 

The  examples  in  Appendix  IV  illustrate  the  computer  program  and 
show  how  it  is  used.  An  example  is  given  here  to  illustrate  the  con- 
cepts of  deterioration  and  reliability.  A "type  1"  or  gravel-surfaced 
facility  is  selected  having  a surface  CBR  greater  than  the  subgrade 
CBR.  The  values  of  the  variables  are  arbitrarily  selected  for  illus- 
tration and  are  shown  in  Table  5.  In  this  example,  the  object  is  to 
evaluate  the  rut  depth  and  reliability  at  all  repetition  levels  up  to 
10,000. 

The  allowable  rut  depth,  RDA  , selected  is  2 in.  Eleven  iterations 
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TABLE  5-— VARIABLES  USED  IN  EXAMPLE  PROBLEM 


Index 

Variables 

Mean 

Variance 

1* 

P 

15  kips 

1423,333 

2 

*P 

70  psi 

58.33 

3 

t 

9 in. 

1.21 

1+ 

C1 

25 

6 

5 

C2 

10 

2 

6a** 

R1 

100 

10 

6b 

R2 

1,000 

100 

6c 

R3 

2,000 

200 

6d 

Rl* 

3,000 

300 

6e 

R5 

14,000 

1400 

6f 

R6 

5,000 

500 

6g 

rt 

6,000 

600 

6h 

R8 

7,000 

700 

6i 

R9 

8,000 

800 

6J 

R10 

9,000 

900 

6k 

R11 

10,000 

1,000 

* Values  of  variables  1-5  used  for  all  iterations. 

**  Value  of  variable  6 changed  for  each  iteration. 

of  the  problem  are  performed,  each  time  changing  the  value  of  the  num- 
ber of  repetitions,  R . The  results  are  portrayed  graphically  in 
Fig.  56,  where  the  values  of  rut  depth  in  inches  and  the  reliability  are 
shown  as  a continuous  function  of  the  number  of  15-kip  repetitions. 


91 


The  information  shown  in  Fig.  56  gives  the  designer  or  evaluator 
an  illustration  of  the  effect  of  repetitions  upon  both  rut  depth  and 
reliability.  Although  the  maximum  allowable  rut  depth  is  not  exceeded, 
it  has  an  associated  r iability  of  only  0.5.  If  a higher  reliability 
level  were  required,  structural  change  in  the  facility  would  be  neces- 
sary. In  this  example,  an  increase  in  the  thickness,  t , and/or  the 
CBR,  , of  the  gravel-surfaced  course  would  bring  about  the  desired 
change.  Additional  iterations  of  the  problem  using  new  values  of  t 
and  would  be  required  for  evaluation.  The  data  generated  would 

provide  for  additional  rut  depth  and  reliability  relations  to  be 
plotted  on  Fig.  56. 

Although  this  example  problem  illustrates  how  the  rate  of  deteri- 
oration and  the  reliability  can  be  evaluated,  it  also  illustrated  the 
overall  potential  of  the  DAS,  as  follows: 

1.  A quantitative  time  and  use  rate  of  change  is  developed  such 
that  the  engineer  is  no  longer  constrained  by  failure  point  design 
criteria. 

2.  Limiting  serviceability  criteria  can  be  selected  that  best 
suit  an  engineering  requirement. 

3.  Any  level  of  reliability,  or  degree  of  conservatism,  can 
be  selected  to  meet  the  needs  of  the  designer. 

U . Quantitative  bases  for  programming  of  maintenance  and  repair 
exist  as  a result  of  the  deterioration  analysis  capability. 

5.  Quantitative  differential  analysis  is  possible  because  the  de- 
gree of  deterioration  induced  by  different  quantities  and  types  of  ve- 
hicles can  be  determined. 
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FIG.  56. —EXPECTED  VALUE  AND  VARIANCE  OF  RUT  DEPTH,  AND  ASSOCIATED 
RELIABILITY  VERSUS  NUMBER  OF  15-KIP  ESWL  REPETITIONS 
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6.  Military  planning  is  enhanced  through  the  ability  to  quantify 
the  state  of  deterioration  at  any  anti'  pated  level  of  traffic  usage. 

7-  Thickness  and  strength  design  of  facilities  can  be  optimized 
with  respect  to  degree  of  deterioration,  reliability,  and  use  of  avail- 
able construction  materials. 

1 
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CONCLUSIONS  AND  RECOMMENDATIONS 
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Conclusions 

The  conclusions  drawn  as  a result  of  this  study  are  as  follows: 

1.  The  hypothesis  that  effective  pavement  life-cycle  management  can 
he  achieved  through  the  use  of  deterioration  and  reliability  concepts  has 
been  investigated  and  proven. 

2.  Models  were  developed  that  effectively  portray  the  deterioration 
of  a facility  and  assess  its  reliability  in  terms  of  the  rutting  mode  of 
det  er iorat ion . 

3.  The  DAS  is  a first-generation  system  that  can  be  considered 
validated  for  use  to  the  extent  of  the  current  CE  design  procedure. 

U.  The  deterioration  and  reliability  models  show  high  correlation 
and  small  residual  error  and,  therefore,  when  combined  to  form  the  heart 
of  the  DAS,  should  provide  for  effective  rut  depth  prediction  and  relia- 
bility assessment. 

5.  The  DAS  can  be  used  for  design  and  evaluation  , optimization, 
differential  analysis,  planning,  and  military  operations  as  described 
in  the  section  entitled  "Utilization  of  the  DAS." 

6.  The  DAS,  as  a first-generation  system,  provides  a basis  for 
the  development  of  a complete  life-cycle  management  system  for  all  modes 
of  deterioration  pertinent  to  all  pavement  types  through  expansion  and 
validation  as  data  are  made  available. 

7.  The  DAS  can  be  used  in  its  present  form  for  relative  differ- 
ential analysis  on  roads  where  damage  incurred  by  various  vehicle  types 
must  be  determined  as  a basis  for  cost  assessment. 
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Recommendat ions 


The  scope  of  this  study  and  the  results  obtained  are  limited  in 
terms  of  intended  use.  Although  the  results  can  be  directly  utilized  for 
some  purposes,  the  basic  intent  is  to  establish  the  fact  that  deteriora- 
tion and  reliability  concepts  can  play  a vital  role  in  improving  the 
state  of  the  art  in  pavement  design,  evaluation,  maintenance,  and  overall 
life-cycle  management.  Additional  validation  and  expansion  are  required 
for  verification  and  improvement  of  the  overall  system. 

The  following  specific  recommendations  are  considered  appropriate: 

1.  The  DAS  should  be  validated  for  use  in  terms  that  satisfy  the 
needs  of  the  using  agency. 

2.  The  DAS  should  be  employed  on  a trial  basis  by  the  CE,  FS,  and 
other  appropriate  agencies  to  increase  awareness  and  determine  possible 
benefits  that  can  be  derived. 

3.  The  current  field  evaluation  programs  should  be  continued  and 
improved  to  provide  more  closely  controlled  data  and  expanded  to  in- 
corporate all  of  the  more  important  deterioration  modes. 

!+.  The  present  research  programs  should  be  continued  and  expanded 
in  scope  and  level  of  intensity  to  effectively  provide  for  development 
of  a fully  validated  and  comprehensive  system  within  this  decade.  A 
significant  portion  of  these  expanded  programs  should  include  the 
investigation  of  the  effects  of  vehicle  dynamics  as  well  as  the  investi- 
gation of  the  correlation  among  steady  state  tangent  operations,  curve 
operations,  and  acceleration-deceleration  operations. 
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6.2 

35 

0.15 

10000 

100 

1 1 

132 

6.2 

141 

0.20 

1 0000 

100 

1 1 

132 

6.2 

353 

0.20 

10000 

100 

1 1 

132 

6.2 

706 

0.19 

10000 

100 

1 1 

132. 

6.2 

1765 

0.20 

1 0000 

100 

1 1 

132 

6.2 

3530 

0.30 

10000 

100 

1 1 

132 

6.2 

6001 
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F5VL 

TP 

TH1 

TH2 

CBR1 

CBP2 

REP 

.20 

20000 

50 

3 

22 

54 

6.5 

254 

.90 

20000 

50 

3 

22 

54 

6.5 

1247 

1 .10 

20000 

50 

3 

22 

54 

6.5 

8951 

1 .20 

20000 

50 

3 

22 

54 

6.5 

12700 

.20 

20000 

50 

3 

17 

53 

5.5 

254 

.80 

20000 

50 

3 

17 

53 

5.5 

1247 

1 .30 

20000 

50 

3 

17 

53 

5.5 

8951 

1 .60 

20000 

50 

3 

17 

53 

5.5 

12700 

.40 

20000 

50 

3 

11 

51 

4 

254 

1 .00 

20000 

50 

3 

1 1 

51 

4 

790 

1 .60 

20000 

50 

3 

1 1 

51 

4 

1247 

.72 

77400 

100 

3 

21 

84 

4.6 

85 

.96 

77400 

100 

3 

21 

84 

4.6 

170 

I .68 

77400 

100 

3 

21 

84 

4.6 

1222 

2 . : 6 

77400 

100 

3 

21 

84 

4.6 

4080 

1 .68 

07000 

100 

3 

21 

73 

3.9 

1036 

2.52 

87000 

180 

3 

21 

73 

3.9 

1647 

.4? 

75000 

290 

3 

21 

100 

4.2 

137 

1 .44 

75000 

290 

3 

21 

100 

4.2 

328 

.50 

200000 

150 

9.5 

55.5 

65 

10 

10 

1 .00 

200000 

150 

9.5 

55.5 

65 

10 

50 

1 .55 

200000 

150 

9.5 

55.5 

65 

10 

210 

2.10 

200000 

150 

9.5 

55.5 

65 

10 

530 

2.50 

200000 

150 

9.5 

55.5 

65 

10 

2100 

.40 

200000 

150 

16 

50 

65 

10 

20 

I .00 

200000 

150 

16 

50 

65 

10 

110 

2 .40 

200000 

150 

16 

50 

65 

10 

1060 

2.40 

200000 

150 

16 

50 

65 

10 

2100 

.43 

64000 

200 

1 

14 

100 

47 

628 

.52 

64000 

200 

1 

14 

100 

47 

1727 

.83 

64000 

POP 

1 

14 

100 

47 

3925 

1 .04 

64000 

200 

1 

14 

100 

47 

6280 

1 .16 

64000 

200 

1 

14 

100 

47 

9420 

1.19 

64000 

200 

1 

14 

100 

47 

14353 

.14 

15000 

45 

2 

10 

35 

5 

24 

.24 

15000 

45 

2 

10 

35 

5 

48 

.46 

15000 

45 

2 

10 

35 

5 

299 

.55 

15000 

45 

2 

to 

35 

5 

538 

.81 

15000 

45 

2 

10 

35 

5 

777 

.91 

15000 

45 

2 

10 

35 

5 

892 

.40 

Q4  000 

190 

3 

25 

100 

5.4 

74 

.60 

94000 

190 

3 

25 

100 

5.4 

592 

.80 

94000 

190 

3 

25 

100 

5.4 

2553 

1 .30 

94000 

190 

3 

25 

100 

5.4 

4625 

1 .60 

94000 

190 

3 

2.5 

too 

5.4 

6771 

.30 

94000 

190 

3 

25 

100 

5.4 

74 

.70 

94000 

190 

3 

25 

100 

5.4 

592 

1 .20 

94000 

190 

3 

25 

100 

5.4 

6771 

.30 

94000 

190 

3 

25 

100 

3.8 

74 

.60 

94000 

190 

3 

25 

100 

3.8 

592 

1 .00 

94000 

190 

3 

25 

100 

3.8 

2553 

1 .20 

9400O 

190 

3 

25 

100 

3.8 

4625 

1 .30 

94000 

190 

3 

25 

100 

3 .8 

6771 

1 .40 

94000 

190 

3 

25 

too 

3.8 

14467 

119 


TWO-LAYER  FLEXIBLE  PAVEMENT  DATA 


TV0L4Y 

corn 

.38 

94000 

.80 

94000 

1.10 

94  000 

1 .40 

94000 

1 .v',0 

94000 

1 .90 

94000 

.4*; 

94000 

.70 

94000 

I .00 

94000 

1 .40 

94000 

1.40 

94  000 

1 .80 

94000 

.40 

94000 

.80 

94000 

1 .58 

94000 

.30 

94000 

.80 

94  000 

1 .30 

94000 

.30 

94000 

.60 

94000 

1 .10 

94000 

.40 

1 13000 

.40 

1 13000 

.70 

113000 

.40 

1 13000 

.70 

1 13000 

1 .20 

1 13000 

,5P 

1 13000 

.80 

1 13000 

I .40 

1 13000 

.40 

1 13000 

.60 

1 13000 

1 .10 

1 13000 

1 .90 

1 13000 

.40 

1 13000 

.70 

1 13000 

1.30 

1 13000 

2.50 

1 13000 

.40 

1 13000 

.88 

1 13000 

1 .20 

I 13000 

1 .90 

113000 

.90 

1 13000 

1 .70 

1 13000 

.80 

1 13000 

1 .60 

1 13000 

.80 

113000 

.15 

15000 

.20 

15000 

.25 

15000 

.30 

15000 

.40 

15000 

.50 

15000 

.55 

15000 

.60 

15000 
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190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

190 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

258 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250. 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

250 

3 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 
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25 

100 

3.8 

25 

108 

3.8 

25 

100 

3.8 

25 

1P0 

3.8 

25 

100 

3 .8 

25 

100 

3.8 

25 

180 

3.8 

25 

100 

3.8 

25 

100 

3.8 

25 

100 

3.8 

25 

100 

3.8 

25 

100 

3.8 

25 

96 

4.9 

25 

96(® 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

96 

4.9 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

100 

4 

25 

108 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

3.2 

25 

100 

5.2 

25 

100 

5.2 

25 

100 

5.2 

25 

100 

5.2 

25 

100 

5.2 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 

6.5 

55 

35 
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74 

592 

2553 

4625 

6771 

14467 

74 

592 

2553 

4625 

6771 

14467 

74 

592 

2553 

74 

592 

2553 

74 

592 

2553 

74 

222 

592 

74 

222 

592 

74 

222 

592 

74 

222 

592 

1147 

74 

222 

592 

1147 

74 

222 

592 

1147 

74 

222 

74 

222 

74 

22 

42 

100 

190 

380 

960 

1910 

3820 
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.65 

15000 

.15 

15000 

.20 

15000 

.35 

15000 

.50 

I5O0O 

.60 

15000 

.85 

15000 

1.05 

15000 

1 .25 

15000 

1 .40 

15000 

.15 

15000 

.25 

15000 

.43 

15000 

.55 

15000 

.75 

15000 

1 .00 

15000 

1.25 

15000 

1 .50 

15000 

1.55 

I5000 

.15 

15000 

.20 

15000 

.20 

I5O0O 

.25 

15000 

.25 

15000 

.35 

15000 

.40 

15000 

.45 

15000 

.45 

15000 

.10 

15000 

.10 

15000 

.15 

I5000 

.20 

150O0 

.25 

15000 

.50 

15000 

.55 

15O0O 

.70 

15000 

.85 

150O0 

.35 

I50P0 

.45 

15000 

.50 

15000 

.65 

15000 

.95 

15000 

1 .20 

15000 

I .50 

15000 

1 .80 

15OO0 

2.00 

15000 

.20 

37000 

.45 

37000 

.60 

37000 

.95 

37000 

1 .20 

37000 

1.55 

37O0P 

1 .90 

370OO 

2.25 

37000 

2.50 

37000 
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45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 
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45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 
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45 

3 

45 

3 

45 

3 

45 

3 

45 

3 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

45 

1.5 

04/21 

/78 

6.5 

55 

35 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

6.5 

53 

29 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

59 

27 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

55 

35 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

5 

77 

38 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

6.5 

61 

44 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 

8.5 

60 

34 
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5690 

20 

40 

100 

190 

380 

960 

1910 

3820 

6690 

20 

40 

too 

190 

380 

960 

1910 

382? 

6690 

20 

40 

100 

190 

380 

960 

1910 

3820 

6690 

20 

40 

100 

190 

380 

960 

1910 

3820 

6690 

20 

40 

100 

190 

380 

960 

1910 

3820 

6690 

15 

30 

80 

160 

320 

800 

1590 

3180 

5570 
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.25 

37000 

45 

3 

7 

88 

32 

15 

.40 

37000 

45 

3 

7 

88 

32 

30 

37000 

45 

3 

7 

88 

32 

80 

.90 

37000 

45 

3 

7 

88 

32 

160 

1 .05 

37000 

45 

3 

7 

88 

32 

320 

1 .50 

37000 

45 

3 

7 

88 

32 

800 

1 .80 

37000 

45 

3 

7 

88 

32 

1590 

2.10 

37000 

45 

3 

7 

88 

32 

3180 

2.45 

37000 

45 

3 

7 

88 

32 

5570 

.50 

37000 

45 

1.5 

8.5 

63 

44 

15 

.05 

37000 

45 

1.5 

8.5 

63 

44 

30 

1.45 

37000 

45 

1.5 

8.5 

63 

44 

80 

1 .95 

37000 

45 

1.5 

8.5 

63 

44 

160 

2.50 

37000 

45 

1.5 

8.5 

63 

44 

320 

3.00 

37000 

45 

1.5 

8.5 

63 

44 

800 

3.50 

37000 

45 

1.5 

8.5 

63 

44 

1590 

4.00 

37000 

45 

IN5 

8.5 

63 

44 

3180 

4.15 

37000 

45 

1.5 

8.5 

63 

44 

5570 

.20 

37000 

45 

3 

7 

60 

34 

15 

.35 

37000 

45 

3 

7 

60 

34 

30 

.50 

37000 

45 

3 

7 

60 

34 

80 

.75 

37000 

45 

3 

7 

60 

34 

160 

1 .00 

37000 

45 

3 

7 

60 

34 

320 

1 .45 

37000 

45 

3 

7 

60 

34 

800 

1 .60 

37000 

45 

3 

7 

60 

34 

1590 

2 .00 

37000 

45 

3 

7 

60 

34 

3180 

2.10 

37000 

45 

3 

7 

60 

34 

5570 

.65 

37000 

45 

1.5 

8.5 

41 

29 

30 

1 .20 

37000 

45 

1.5 

8.5 

41 

29 

80 

1 .05 

37000 

45 

1.5 

8.5 

41 

29 

160 

2.15 

37000 

45 

1.5 

8.5 

4 1 

29 

320 

2 .P5 

37000 

45 

1.5 

8.5 

41 

29 

800 

3.15 

37000 

45 

1.5 

8.5 

41 

29 

1590 

3.55 

37000 

45 

1.5 

8.5 

41 

29 

3180 

3.90 

37000 

45 

1.5 

8.5 

41 

29 

5570 

.40 

37000 

45 

3 

7 

48 

27 

8 

.90 

37000 

45 

3 

7 

48 

27 

15 

1 .3  5 

37000 

45 

3 

7 

48 

27 

30 

2.10 

37000 

45 

3 

7 

48 

27 

80 

2.00 

37000 

45 

3 

7 

48 

27 

160 

3.35 

37000 

45 

3 

7 

48 

27 

320 

4.00 

37000 

45 

3 

7 

48 

27 

800 

4.50 

3 7000 

45 

3 

7 

48 

27 

1590 

4.95 

37000 

45 

3 

7 

48 

27 

3180 

5.05 

37000 

45 

3 

7 

48 

27 

5570 

.25 

50000 

50 

1.5 

8.5 

65 

34 

15 

.40 

50000 

50 

1.5 

8.5 

65 

34 

30 

.65 

50000 

50 

1.5 

8.5 

65 

34 

R0 

1 .00 

50000 

50 

1.5 

8.5 

65 

34 

160 

1 .30 

50000 

50 

1.5 

8.5 

65 

34 

320 
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50000 

50 

1 .5 

8.5 

65 

34 

800 

2.00 

50000 

50 

1.5 

8.5 

65 

34 
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2.45 
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1 

12 

100 

29 

92 

13 

12560 

.48 

85000 

200 

1 

12 

100 

29 

92 

13 

14353 

.06 

75000 

200 

1 

12 

100 

15 

63 

12 

63 

.29 

75000 

200 

l 

12 

100 

15 

63 

12 

628 

.56 

75000 

200 

1 

12 

100 

15 

63 

12 

6280 

.59 

75000 

200 

1 

12 

100 

15 

63 

12 

14353 
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APPENDIX  IV. —INPUT,  OUTPUT,  AND  PROGRAM  LISTINGS 


The  RUTDEP  program  was  written  in  Fortran  IV  computer  language  and 
was  run  using  the  WES  600  time-sharing  computer  system.  A complete  list- 
ing of  the  program  is  included.  With  minor  modifications,  the  RUTDEP 
program  should  he  adaptable  to  any  computer  system  if  the  following 
conditions  are  satisfied: 

1.  A time-sharing  computer  system  that  will  allow  Fortran  IV 
computer  code  must  be  used. 

2.  A remote  terminal  is  necessary  to  access  the  time-sharing 
computer  system. 

3.  The  RUTDEP  program  must  be  stored  on  disc  in  the  time-sharing 
system. 

4.  The  ability  to  create  and  store  input  data  files  on  disc  in  the 
time-sharing  system. 

When  these  conditions  are  satisfied,  the  program  can  be  executed  as 
shown  in  the  following  example.  For  ease  in  running  the  RUTDEP  program, 
the  input  data  are  entered  in  a step-by-step  process  in  response  to 
questions  printed  by  the  remote  terminal.  Included  with  the  example 
problem  shown  below  is  a discussion  of  each  step.  The  responses  of  the 
computer  user  are  underlined. 

The  initial  step  is  to  access  the  computer  system  and  call  up  the 
RUTDEP  program.  With  the  WES  600  computer  system  this  is  done  as 
follows : 

***  WES-TSS  NOTIFIED  *** 

HIS  SERIES  600  ON  01/31/78  AT  10.750  CHANNEL  6575 
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APPENDIX  IV. — INPUT,  OUTPUT,  AND  PROGRAM  LISTINGS 


USER  ID  -R0SF111 
PASSWORD— 

XXXXXXXXXXX 

SYSTEM  ? FORTRAN 

OLD  OR  NEW-OLD  RUTDEP 

READY 

*RUN 

After  the  RUTDEP  program  is  called,  the  run  command  is  given  and 


the  computer  asks  for  the  type  of  problem  to  be  run. 

TYPE  PROBLEM: 

UNSURFACED  C1>C2  = 1 
UNSURFACED  C1<C2  = 2 
ASPHALT  WO/SUBBASE, 8 VAR.  = 3 
ASPHALT  W/SUBBASE,10  VAR.  = U 


=k 

A listing  of  the  four  types  of  problems  is  printed  out  with  an  iden- 
tification number  (l,  2,  3,  M for  each  type.  The  appropriate  number  is 
then  entered  for  the  type  of  problem  the  user  requires.  In  this  example 
an  asphalt  pavement  with  a subbase  was  required,  so  the  number  "U"  was 
entered. 

INPUT  - ALLOWABLE  RUT  DEPTR. 

= 1.0 


The  next  item  of  input  is  the  allowable  rut  depth,  which  is  used  in 
computing  the  reliability  statistic.  An  allowable  rut  depth  of  1.0  inch 
was  selected. 


INPUT  MODE  - FILE(l), KEYBOARDS), COMBINATION^)? 
= 2 


The  next  question  asked  by  the  computer  is  how  the  actual  problem 
data  will  be  entered  into  the  computer;  by  a data  file,  by  the  keyboard 
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on  the  terminal,  or  by  a combination  of  these  two.  The  computer  user 
selects  one  of  the  three  modes  by  entering  the  identification  number 
(1,  2,  3). 

Depending  upon  the  mode  selected  the  computer  asks  one  of  two  ques- 
tions. If  all  of  the  data  are  to  be  entered  by  the  file  mode  (l)  or  if 
part  of  the  data  are  to  be  entered  by  a file  (3),  the  computer  responds 
by  asking  for  the  name  of  the  data  file.  Instructions  on  how  to  create 
an  input  data  file  will  be  discussed  later.  In  this  example,  the  second 
input  mode  was  selected,  which  is  input  by  the  keyboard.  Since  no  data 
file  is  required  with  this  mode,  the  computer  did  not  ask  for  a file  name 
but  went  to  the  next  step,  asked  for  the  input  stat  code  of  the  vari- 
ables, and  listed  the  variables  in  order. 

INPUT  - STAT  CODE  - COMPUTED(O),  GIVEN(l) 

P TP  R T1  T2  Cl  T3  C2  C3 
=111111111 

A zero  or  one  should  be  entered  on  the  keyboard  by  the  user  for 

each  variable  listed.  The  number  zero  instructs  the  computer  to  compute 

the  mean  and  variance  of  the  variable  from  the  data  file,  while  the 

number  one  tells  the  computer  that  the  mean  and  variance  will  be  typed 

in  from  the  terminal.  Since  no  data  file  was  used  for  this  example 

problem,  the  number  one  was  entered  for  all  nine  variables.  The  computer 

then  asks  for  the  mean  and  variance  of  the  first  variable.  The  mean  and 

variance  may  be  separated  by  a comma  or  a space. 

INPUT  MEAN  & VARIANCE  FOR  - P 
=18000.  500000. 

When  the  data  are  entered  for  the  ESWL  (P)  the  computer  asks  for 
the  input  data  of  the  second  variable. 
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I 


INPUT  MEAN  & VARIANCE  FOR  - TP 
=80.  60. 


After  the  user  responds  to  this  question,  the  computer  asks  for  the 
mean  and  variance  of  the  remaining  variables  one  at  a time. 


INPUT  MEAN 

& 

VARIANCE 

FOR  - 

R 

=40000.  1000000. 

INPUT  MEAN 

& 

VARIANCE 

FOR  - 

T1 

=2.3  0.5 

INPUT  MEAN 
=6. ,1. 

& 

VARIANCE 

FOR  - 

T2 

INPUT  MEAN 

& 

VARIANCE 

FOR  - 

Cl 

=50. ,10. 
INPUT  MEAN 

& 

VARIANCE 

FOR  - 

T3 

=6.  1. 
INPUT  MEAN 

& 

VARIANCE 

FOR  - 

C2 

=25.  L_ 
INPUT  MEAN 

& 

VARIANCE 

FOR  - 

C3 

=10. ,2. 


In  entering  data  by  either  a data  file  or  by  the  keyboard,  the  ESWL 
(P)  should  always  be  in  pounds,  the  tire  pressure  (TP)  should  be  in  psi, 
and  all  thicknesses  should  be  in  inches.  After  the  mean  and  variance 
for  the  last  variable  is  entered,  the  computer  executes  the  program  and 
lists  the  output. 

STAT 


VARIABLE 

MEAN 

VARIANCE 

CODE 

P 

18000.00 

500000.00 

1 

TP 

80.00 

60.00 

1 

R 

40000.00 

1000000.00 

1 

T1 

2.50 

0.50 

1 

T2 

6.00 

1.00 

1 

Cl 

50.00 

10.00 

1 

T3 

6.00 

1.00 

1 

C2 

25.00 

7.00 

1 

C3 

10.00 

2.00 

1 

The 
the  mean 


first  item  of  input  is  a listing  of 
and  variance  of  each  variable.  The 


the  variables  along  with 
stat  code  listing  of  the 
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output  indicates  if  the  mean  and  variance  of  each  variable  was  computed 
(0)  or  given  by  the  user  (l). 


MEAN  VALUE  OF  RUT  DEPTH  = 0. 865862 
EXPECTED  VALUE  OF  RUT  DEPTH  = 0.  90011+6 
EXPECTED  VARIANCE  OF  RUT  DEPTH  = 0.019569 
RELIABILITY  STATISTIC  = 0.713800 


The  remaining  output  completes  the  problem.  The  computer  then  re- 
peats the  type  of  problem  question. 

TYPE  PROBLEM: 

UNSURFACED  C1>C2  = 1 
UNSURFACED  C1<C2  = 2 
ASPHALT  WO/SUBBASE, 8 VAR.  = 3 
ASPHALT  W/SUBBASE.10  VAR.  = k 


If  another  problem  is  to  be  run,  the  user  enters  the  proper  identifi- 
cation number  as  before.  If  no  more  problems  are  to  be  run,  the  number 
9 should  be  entered  and  the  computer  will  terminate  the  program. 

In  the  above  example,  the  mean  and  variance  data  were  input  from 
the  terminal  keyboard.  If  the  computer  user  had  desired  to  respond  to 
the  input  mode  question  by  selecting  a file  input  mode  (l)  or  a combina- 
tion of  file  and  keyboard  (3),  then  a data  file  must  be  created  and 
stored  before  the  RUTDEP  program  is  run.  The  mean  and  variance  is  then 
computed  from  the  data  points  in  the  data  file.  With  the  file  mode,  the 
data  file  should  be  created  with  the  following  format: 


Line  No. 
Line  No. 
Line  No. 
Line  No. 


N1 

Value  ( 1 ) 
N2 

Value  (2) 


Value  ( 2 ) 
Value  (2) 


Value  (N  ) 
Value  (N2) 


13U 


After  the  data  file  is  created,  it  should  be  stored  under  a unique 
data  file  name,  which  may  have  a maximum  of  six  alphanumeric  character- 
istics. An  example  of  a data  file  for  pavement  type  1 is  shown  below. 

It  was  stored  under  the  file  name  "DATA  1." 

*LIST  DATA1 

100  6 

101  lUOOO.  lUUOO.  15000.  15300.  15300.  16000. 

110  6 

111  60.  60.  70.  75.  75.  8o. 

120  5 

121  U200.  1*800.  5000.  5500.  5500. 

130  7 

131  7-5  8.0  8.5  9.0  9.0  10.0  11.0 

ll*0  1* 

ll+l  21.  25.  27.  27. 

150  3 

151  9.  9.  12. 

Data  points  for  all  six  dependent  variables  were  listed  with  three 
seven  points  per  variable. 

When  test  data  are  available  for  only  some  of  the  variables,  then  a 
combination  of  file  and  keyboard  input  (input  mode  3)  would  be  required. 
The  data  file  would  be  created  in  a format  similar  to  that  of  input 
mode  1.  For  variables  without  data  points,  a one  would  be  entered  for 
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IL  , and  a zero  would  be  entered  on  the  following  line  number  where  the 
values  of  the  data  points  would  normally  be  listed.  An  example  of  this 
type  of  data  format,  which  was  stored  under  the  data  file  name  "DATA  2," 
is  shown  below. 

*LIST  DATA2 

100  6 

101  1^+000.  IMOO.  15000.  15300.  15300.  16000. 

200  1 

201  0. 

300  1 

301  0. 

1*00  7 

1+01  7-5  8.  8.5  9-  9.  10.  11. 

500  1+ 

501  21.  25.  27.  27. 

600  1 

601  0. 


Data  points  for  the  first,  fourth,  and  fifth  variables  are  listed, 
while  a one  and  a zero  are  listed  for  the  other  three  variables.  The 
sequence  of  line  numbers  for  file  "DATA  2"  is  different  from  file 
"DATA  1."  The  sequence  could  have  been  identical  if  desired.  The  only 
rule  that  must  be  followed  is  that  the  line  numbers  must  be  in  ascending 
order. 

In  the  following  examples,  the  three  different  types  of  input  modes 
were  used  to  work  one  problem.  When  the  file  mode  input  (l)  and  combina- 
tion mode  input  (3)  are  selected,  the  computer  asks  for  the  data  file 
name.  In  these  examples,  DATA  1 was  given  for  mode  1,  and  DATA  2 was 
given  for  mode  3.  The  next  question  for  all  three  modes  is  the  stat 
code  for  all  three  modes.  For  mode  1,  a zero  is  entered  for  all  varia- 
bles, and  for  mode  2,  a one  is  entered  for  all  variables.  For  mode  3,  a 
zero  is  given  for  the  variables  with  data  points,  and  a one  is  given  for 
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the  variables  where  the  user  has  to  type  in  the  mean  and  variance. 

Following  the  last  example,  the  computer  again  asked  for  the  type 
of  problem.  A nine  was  entered  and  the  program  terminated. 


HIS  SERIFS  600  ON  04/?  1/7*  AT  10.761 


CHANNEL  5416 


USER  ID  - 

ILLEGAL  ID-RETYPE--P0SF1 1 I 
PASSWOPD-- 

SYSTEM  7F0PT 

ALL  USERS  SFE  INFO  7PEST0PF  FOP  PESTOPES  FFOM  7TPK  556BPI  FIl.SYS  SAVF  TAB 

OLD  OP  NFV-NFV 
READY 

♦OLD  WFSl IB/PAGFF.PeE 

YOU  HAVE  FILE  ACCFSSED  WITHOUT  PEAD  PERMISSION 
SYSTFM  7F0PT  0 WESLIB/PAGEO  ,P 
READY 
♦ RUN 


04/?1/7h  10.7OP 

PROGRAM  PAGEP  --  LISTS  FPOM  ANY  PAGE  TO  THE  END  OF  THE  FILE, 
WITH  OP  WITHOUT  SPECIAL  " NEW  PAGE"  CHAPACTFR. 

(PPOGPAM  PAGES  LISTS  SELECTED  PAGES  ONLY.) 


ENTEP  T IF  YOU  WANT  THE  PAGEP  LISTING  ON  YOUP  TIME-SHARING  TERMINAL 

OP  R IF  vpti  WANT  IT  ON  YOtiP  OFFICE'S  PATCH  PRINTER  - 

- T 

ENTEP  CATALOG  DESCRIPTION  (40  CHARACTERS  MAY)  OP  A 
"CARPI AGE  PETUPN"  IF  IT  IS  IN  YOUR  AFT  OP  USER 
MASTER  CATAl  OG  r 
ENTEP  FILE  NAME  r PUTDEP 

ENTEP  PAGE  SI7F  IN  INCHES17  11 

ENTEP  STARTING  VALUE  OF  PAGE  LOCATION  NUMBER  IN 
FILF  TO  PE  LISTED  = I 

ENTEP  STAPTING  VALUE  FOP  PAGE  NUMBER  TO  BE  PRINTED 
AT  THF  BOTTOM  OF  EACH  PAGE  (OMITTED  IF  7FF0)  r 
ENTEP  YOUP  "MEW  PAGF"  INDICATOR  CHARACTER 
OR  A "CARRIAGE  RETURN"  TO  NOT  USE  ONF  = * 

"«w  page"  indicator  t'ST  be  in  first  ?d  columns  of  lime. 
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RUTDEP  10;50s40  04/PI/7P  FILE  PAGE  NO.  1 


1000C***** PAGE  10 

1 0 1 PC*****  ************************************************************  * 
1020C 

1 03 PC  PROGRAM  PUT  - COMPUTES  PUT  DEPTH  FOP  FOUR  POSSIBLE 
10A0C  PAVEMENT  TYPES. 

1050C 

10GUC  CODED  FOP  WES  GPP  TI MESH APE  COMPUTER  SYSTEM. 

107PC 

10R0C  CODED  BYf  H.  P.  AUSTIN 

1090C  E.  ODOM 

110PC  SOILS  A PAVEMENTS  LABOPATOPY 

111  PC  WATEPWA vp  EXPERIMENT  STATION 

1 1  ?.  PC  SEPTEMBER  1977 

1 1 3PC 

1 1 4 PC ********************** ********************** ****** ****** ********** 
1150  DIMENSION  PD  (3PP ) , P (3PP) , TP (3PP I , T 1 (3PP)  ,T2  (3PP)  ,C  U3P0) 

1 1 GP  DIMENSION  T3 (3PP) , C? (3PP) ,C3 f 3PP) ,F ( 3PP) 

1170  DIMENSION  CASE ( A ) fM ( 1 P> , V( I PI .COMPETE ( IP) 

HEP  DIMENSION  SKIP  (3PP) 

1 1 op  I NTEGF.P  TYPE  COMPUTE  .CASE  .START  ,S  TOP 

12PP  PEAL  'FD.MP.MTP.MTI .MT2.MCI , MT3 ,MC? ,MC3 ,MP 

12  10  FFAl.  M.K 

1220  CHAPACTFP  FNAME*  G,FNAME2*P 

1230  CHAPACTFP  CASE  1*2 ( G) ,C ASE2*2 fG) ,C ASE3*? ( 7) ,C ASFa*? (9 I 

12AP  DATA  CASE/G  G.7,9/  , . 

1250  DATA  CASF 1 / P ","TP",'P  ",“T1  , Cl  , C2  / 

12GP  DATA  CASF2/-P  ^.^TP^.^P  * .^T l' ,”C I ' . 

1270  DATA  CASE3/  P .,.TP.,.P  . . J I. . . T?„ . „C 1 . . „C ‘3  / 

1280  DATA  CASF 4/  P , TP  , P , TI  , T2  , Cl  , T3.C2.C3/ 

129PC 

I3PPC 

1310  I rP 

1320  I01P  CONTINUE 

1 33  PC 

134PC  PEAD  INPUT  DATA 

1350  GO  TO  2P 

13GP  1P2P  CONTINUE 

13  7PC 

13 RPC  PRINT  MEAN  AND  VARIANCE  PESUl TS 

13«P  GO  TO  30 

1400  1030  CONTINUE 

14I0C 

1 42PC  SELECT  AND  SOLVE  COPFFCT  PUT  DEPTH  EPUATION  AND  OUTPUT  RESULTS 
1430  GO  T0(4l,4?,43,44),TYPE 

1440  1040  CONTINUE 

1450C  LOOP  BACK  THROUGH  FOP  ANOTHER  PROBLEM 

1 4G0  GO  TO  1010 

1470  1090  CONTINUE 

14RPC 

1490C  END  OF  PROGRAM 

1500  PRINT," NOPMAl  TERMINATION" 

1510  STOP 

152  PC 

153PC***** — end  of  page  10 — 
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RUTPFP 


CONT 


10:50:40  04/21/78 


FILE  PARE  NO 


2 


I540C* 
I550C 
1560C 
1570  20 
1580C 
1 590C 
1 600 
1610 
1 620 
1630 
16A0 
1650 
1660 
1670C 
1 6P0 
16?0 
1700 
171PC 
I720C 
1 73  OC 
1740 
1750 
1760C 
1770C 
1700 
17O0C 
IPOO  II 
18!0C 
1820C 
1P30C 
1840 
1850 

I860  I 
1P70 

1800  2 
1890 

1900  13 

1910 
1920C 
1930C 
1940C 
1950C 
I960 
1970 
1980 
1990 
2000C 
20  IOC 
2020 
2030 
2040 
2050 
2060 
2070 
2080 


--PAGE  20--- 

PEAR  INPUT  DATA 
CONTINUE 


I -I  1 

PRINT, "TVPF  PROBLEM:" 

PRINT,"  IJNSURFACEP  C I >C2  = l" 

PRINT,"  UNSURFACEP  C1<C?  = 2" 

PRINT,"  ASPHALT  W0/SUBBASE.8  VAR.  : 3 

PRINT,"  ASPHALT  W/SUBPASE.IO  VAR.  = 4 

PFAP (5,P00)TYPE 

IF(TyPE.GT.4)R0  TO  1090 

PRINT, "INPUT  - Al  LOVABLE  PUT  PFPTH." 

PFAP(5,PO0)APP 


PRINT  "INPUT  MOPF  - FILE( 1 ) .KEYBOARD (2) .COMBINATION (3) ?" 
REAP (5,P00)I0 


GO  TO  ( I 1 , 1 3 , I I ) , I 0 
PATA  FILE  Will  PE  USED 
CONTINUE 

PATA  FILE  WHERE  PAVEMENT  DATA  IS  STORED 

PRINT  "DATA  FILE  NAME?" 

READ (5,1) FRAME 
FORMAT (A  6) 

ENCODE ( FNAME2  L2 )"/", FRAME ," t" 

F0PMAT(AI,A6,A  I) 

CALL  ATT ACH  ( 1 0, FN AME2 ,3,0,,) 

CONTINUE 

NVAR  =CASE (TYPE) 

TYPE  A 1 IF  KfAN  A VARIANCE  WILL  BE  FURNISHED  FROM  KEYBOARD 
TYPE  A 0 IF  f AN  A VARIANCE  COMPUTED  IN  SUBROUTINE  STAT. 

IF(TYPF.LT.3)PPINT  821.CASE1 
I F (TYPE ,EP.3)PPINT  823.CASE3 
IF(TYPE.E0.4)PPINT  822.CASE4 
PE AD (5,800) (COMPUTF(II),I I =1 , NVAR) 


DO  199  l - I , NVAR 
1.2  =L 

IF (COMPUTF ( 12  ) ,FP  ,0)CALL  STAT(M,V,L2) 
IF(I0.LT.3)G0  TO  199 
IF (COMPUTE (IP  ) ,E0 ,0)G0  TO  199 
PFAP (10, 800)LINENO. NONE 
ST  APT  z 1 


RUTDFP  CONT 


10;50:40  04/21/78  FILE  PAGF  NO.  3 


2090 
2100  25 

21  IP 
2120  22 
2130 

21 40  21 
2150 
2 1 60 
?!  70 
2180  199 
2190  800 
2200  821 

22  1 03 
2220  822 
22303 
22*0  823 
22503 

22  60C 

22  70C 
2280 
2290 
2300 
2310 
2320 

2330  824 
2340 
2350  201 

23  60 
23  700 


ST0P=8 

CONTI NUF 

IE(N0NE-ST0P)?2,22,2 1 

READ(lP,800)LINEtlO,(SKIP(L3),L3=START,NONF) 
GO  TO  199 

PEAP(10,p00)I.INENO.  (SKIP  (1.3), L3=ST  APT  .STOP) 

S7APT:ST0P-M 

STOPrSTOP+8 

RO  TO  25 

CONTINUF 
FORMAT  (V) 

FOPMAT  ("  INPUT  - STAT  COLE  - COMPUTER  ( 0)  .GIVFN  ( 1 ) " . 
/6(A2 ,2X) ) 

FORMAT ("INPUT  - STAT  CO"F  - COMPUTED (0) .GIVEN ( I )" . 
/9(A2,2X) ) 

FORMAT ("INPUT  - STAT  CODE  - COMPUTED (0)  ,G I VEN ( 1 )"  , 
/7(A2,2X)) 


DO  201  L=  I NVAP 

1F(COMPUT£U.).EO.0)GO  TO  201 
IF(TYPE.LT.3)PPINT  824.CASF  I (1.) 
IF(TvPF.E0.3)PPINT  824tCASE3(L) 
IF(TYPF.E0.4)PPIMT  824.CASE4(L) 

FORMAT ("INPUT  PEAN  3 VAPIANCE  FOP  - " A2) 

PFAD  (5,800)  M(L)  ,V(I.) 

CONTINUF 
RO  TO  1020 

--END  OF  PAGE  20--- 


PUTDEP 


CO  NT 


10:50j40  04/21/78 


FILE  PAGF  NO. 


23  80C  * 

2390C 

2400  30 

24I0C 

2420C 

2430C 

2440 

2450  825 
24803 
2470 
2480 
2490 
2500 
2510 
252  0 
2530C 
2540C 
2550 
2 5600 
2570 
25g0C 
2590C 


826 

204 


PAGE3  0 

CONTINUE 

PRINT  PfAN  3 VAPIANCE  REfUlTS. 

PRINT  825 

FOPMAT(//40X,"STAT*/,"  VARIABLE  l€AN  VAPIANCE 

, 5X, “CODE" ) 

DO  204  Lsl.NVAP 

IF  (TYPF.lt. 3)PPINT  826, CASE  I (L) ,M( L) ,V(L) .COMPUTE (L) 
IF  (TYPE  .EO  .3)  PR  I NT  826.C  ASF3  ( l.)  ,M(1.)  ,V(L)  .COMPUTE  (L) 
IF (TYPE .EP .4 ) PR  I NT  826,CASE4(|.)  ,M(I.)  ,V(L>  ,C  OMPUTE  (L) 
FORMAT  (4X,A2,FI4.2,FI4.2,I8) 

CONTINUE 

RELEASE  INPUT  FILE 
CAll  DETACH (10,,) 

RO  TO  1030 

--END  OF  PAGE  30— 


lUl 


RUTDEP 

CONT  I0j50j40  04/21/7P 

2 6P0C* 

— — m 

PAGE  41--- 

261  PC 

2 62  PC 

Cl 

» C2,  7 VAPIAPl  EG 

2 6?  PC 

2 64  PC 

2650  4J 

COM 

ilNl'E 

266P 

rPztie  I) 

267P 

l'¥>  rMP  / I 000. 

2680 

VP-V(|  ) 

2690 

VP  =VP/| 000000. 

27PP 

(1TP  :C(2  ) 

271P 

VTP  :V  (2  ) 

272P 

-n<  3 > 

2730 

VP  :V  ( 3 ) 

274p 

MTI:M<4) 

2750 

VT  | “V  (A  ) 

276P 

MCI :M<5) 

277P 

VC1 :V(5) 

27g0 

MC2:Mf 6) 

2790 

VC2 :V (6 ) 

2800 

K :P  . 1 74  | 0 

2PIP 

A r0 . 4 7 0 7 

2 P2P 

P :0 .5695 

2830 

C r0.2476 

2S4P 

D:-2 .0020 

2850 

Er-P.9335 

2B6P 

F:-P.28A8 

2g70 

PI  rK*  MP**  A*  MTP**B*  ,f»**C 

2PPP 

P2:  CA1.0G  1 0 ( HI  1 ))**D 

2890 

P3:MCI**1*  Jtf2**F 

2900 

P aP |*P2*P3 

29  IP 

P2:(A*(A-1)*VP)/MP**2 

292  P 

TP2:<P*CP-I  >*VTP  )/P>TP**2 

2930 

P2 : (C*  (C  - 1 „)*  VP)  /!*|R**2 

2940 

T12A  -tf.4343*D*  VT1 

2950 

T 12B  :MT  1**2*  A1.0G  1 0 ( TT  1 ) 

2960 

T12Ct(0.4343*(D-1.))/ALOG10<MT1) 

2970 

TI2:TI2A/TI2B*(T12C-I  .) 

29RP 

CI2:(F*  (E-l  .)*VCI  )/MCI**2 

2990 

C22:<F*CF-I  .>*UC2)/MC2**2 

3000 

EPD:P+(0.5*P*(P2.+  TP2+P2+T12+C  12+C22)  ) 

3010 

P 1 : (A**2*  VP  ) /MP**2 

302P 

TP  1 : (B**2*  VTP ) /MTP**2 

3030 

PI  :(C**2*  W>/MP**2 

3040 

T 1 1 A :<0.4343*D)**2*VT1 

3050 

Tl  IP  :(HTI  )**2*(AL0G  10CMTI  >>**2 

3060 

T1 1 =T 1 1A/T1  IB 

3070 

Cl  I :E*  *2*  WC 1 /MC 1 **2 

3080 

C2 l :( F**2* VC2 ) /MC2**2 

3e9P 

VRD:P**2*  (P  l + TP  l+R  1+  Tl  1+C  1 I+C2  1 )-(0.2 

3I0PA 

CP  2**2+  TP2.**2+P2**2+  Tl  2**2+C  1 2**2 

3110 

WPJTE(6,904)P 

3I2P 

WPITEC6.903)  EPD 

3I3P 

WPITFC6,902) VPD 

3140  902 

FOPMATC IPX, "EXPECTED  VAPIANCE  OF  PUT 

FILE  PAGE  NO. 


DEPTH:", FIG. 6//> 


1U2 


J 


RUTDFP 


CONT 


10:50  :4p  04/21/78 


FIl.F  PAGE  NO 


6 


3150  903 

FORMAT  ( l 0X  , " FXPECTFD  VAl.UF  OF  PUT  DEPTH 

=",F16.6//> 

3160  904 

FOPMAT(//10X,“MEAN  vai.ue  OF  PUT  DEPTH 

s" , FI  6.6//) 

3170 

PF1  = (APD  -F.PD  ) /VPD**0.5 

3180 

WRITE ( 6, 9 1 9 )PFl 

3190  919 

FOPM AT (IPX, "RELIABILITY  STATISTIC 

, FI  6.6//) 

3200C 

3210 

no  to  1040 

322CC 

323  0C 

324  0C 

--END  OF  PACE  4j 

RUTDFP  CONT 


10:50:40  04/21/78  F1LF  PAGE  *10.  7 


3250C* 
32  60C 
32  7 PC 
3280C 

32  9 PC 
3300  J>2 
331P 
3320 
3330 
3340 
3350 

33  60 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 


Cl  « C2,  7 VARIABLES. 


CONTINUE 

MP=M(  1) 

VP  :V  ( I ) 
rfTPrm?) 

VTP  =V(2> 

1R =M(3) 

VP  -V  (3 ) 
f*TT  I =M(  4 ) 

VTI=V(4) 

MCI =M(5> 

VCI =V(5) 

MC2=N(6) 

VC2  = V ( 6) 

K =0.1  1009 
A =0.4925 
B =0.854g 

C =0.5018 
0=0.4293 
E=-l .9773 
F=-l .2015 

P I r K*  MP**A*  MTP**B*  i*F**C 
P2.=  (Al  OC  10  ( Ml  1 ) )**  D 
P3=MC1**  * MC2**F 
P dQ  1*02*03 

P2.=  (A*(A-1)*VP)/MP**2 
TP2=(B* (B-l >*  VTP)/MTP**2 
P2=(C*  (C  - I . >*  W>  /MP**2 
T12A  =0.4J43*  r*VT  I 
TI2B=MTI**2*AlOCI0fMTI  ) 

TI2C -(0.4  34  3*  (0-1  .11  /ALOGl0(MTt  1 
T 12  =T  I2A  /T  12B*  ( T 12C  - 1 . ) 

CI2  = (F*  (F.-l  .)*VCn/MCl**2 
C22  = ( F* ( F- I . )*  VC2 ) /MC2**2 
FFD =P+(0.5*P* (P2+  TP2+P2+  T1 2+C I 2+C22  ) ) 


Ri'TDEP 


CONT 


I0j  50:40 


04/21/78 


FILE  PAGE  NO 


3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730A 

3740 

3750 

3760 

3770 

3780 

3790C 

3800 

38I0C 

3 82 PC 

3830C 


P|:(4**?*VP)/MP**? 

TP  I = (B**2*  VTP  ) /t*TP**2 
PI =(C**2*  VP)  /MP**2 
T I 14  -(0..4  54  3*n)**?*VTl 

ti  ib  -(nrn**?*u>i  OGiPCfn  )>**2 
T 1 1 =T  1 14  /T  1 IB 
Cl  1 =E*  *?*  VC  I /MC I **2 
C?  I : ( F**2*  VC?  ) /PC?**? 

VPD=P**?* (P  1+  TP  1+P  1+  TI  1+C1 I+C? 1)- (0.25*0**?)* 
(P?**?+  TP2**?+P?**?+TI  ?.**2+C  I ?**2+C??**? ) 
WRITF (6,90^ )0 
WPITF (6  903)  O'P 
WPITF(6,90?) VPP 
PEI  = (4PP  -F.FP  ) /VPD**0 .5 
WRITE  C 6t  9 1 9)PFL 

GO  TO  1040 


— -END  OF  P4GF  4? 


RUTDEP  CONT 


10:50:40  04/21/78 


FILE  PAGE  NO.  9 


3g4flC*  PAGE  A3 

3850C 

38SPC  ASPHALTIC  CONCRETE  WITHOUT  SUBGPADE  8 VARIABLES. 

3870C 

388CC 

3890  A3  CONTI  lUF. 

3900  MP  =1*1(1) 

3910  !*P=MP/I000. 

3920  VP  =1/(1) 

3930  VP =VP/ 1000000. 

39A0  riTPrCK?) 

3950  VTPrV(?) 

3960  r*F=M(3) 


3970 
3980 
3990 
A000 
A010 
4020 
A030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4 1?  0 
4130 
4140 
4150 
4160 

41  70 
4180 
4190 
4200 
4210 
4220 

42  3 0 
4240 
4250 
4260 
4270 
4280 

42  9 0 
4300 
4310 
4320 
4330 
4340 
4350 
4360 

43  70 
4380 


VR=V<3) 

’1Tl=f*l(4) 

VTI = V(4) 

1T2  =1*1(5) 

V^r  V(5) 

MCI =M(6) 

VC  1 =V (6) 

MC2=M( 7) 

VC2=V(7) 

K =1 .9431 
A =1  .3127 
B =0.0499 
C =0.3240 
D=-3 .4204 
E=-l .6877 
F=-0. 1 156 

01  =K*flP**A*  f*rTP**B*  ,*P**C 

P2  = (Al  OG  10(  1 .25**TT  1 + MT2  ))**!) 

03  =MC  1**  B*  MC2**  F 
0 1*02*03 

P2=(A*(A-I)*VP)/MP**2 
TP 2= (B* (B -1 )*  VTP ) /l*ITP**2 
P2  = (C*(C-1  ,)*VP)/MP**2 
TI2A  =1  ,25**2*0. 434  3*  D*  VT I 

TI2B  =(  1 .25*  IT  1+  WT2  )**2*Al.0G  10(1 .25*  IT  I + MT2  ) 
TI2C=(  (0.4  343*  (D-l  , ) ) /A LOG  1 0(  I .25*  IT  1 + W2  > ) 
TI2=TI2A/TI2B*(T12C-1 .) 

T22A  =0.4343*D*  VT2 
T22=T22A/TI2B*(TI2C-I .) 

C12  = (E>  (E-l  ,)*VCI  ) /PIC  1**2 

C22  = ( F*  (F-l  ,)*VC2)/MC2**2 

EPD  =0+  (0 .5*0*  (P2+  TP2+P2+  Tl  2+  T22+C  1 2+C22 ) ) 

PI=(A**2*VP)/flP**2 

TP  1 = (B**2*  VTP ) /MTP**2 

PI  = (C**2*V4’)/f1R**2 

Tl  IA  =1 .25**2*  (0.4343*  T))**2*  VTI 

Tl  IB=(1.25*firi+HT2)**2*(ALOGI0(l  .25*  1TI  + MT2  ) )**2 

Tl I =T I IA/TI IB 

T2  IA  =(0.4343*D)**2*VT2 

T2  I =T2  IA/TI  IB 

Cl l=F**2*VC I /MCI**2 

C2 1 = ( F**2* VC2 ) /l*IC2**2 


1U5 


RUTDEP 


10»50»40  P)4/?1/7P 


FI1  E PAGE  NO. 


4390 

4400A 

4410 

4420 

4430 

4440 

4450 

446PC 

4470 

4480C 

4490C 

4500C 


VPP:C**J>*  (P  1+ TP  l+R  1 + T1  I + T2I+C1  1-t-C?  I )-  (0.25*0**?  >* 

(P2**?+  TP  T1  2**2+  T?2**2+C  1 ?**2+C??**2 ) 

WPITECfi, 904)0 
WRITF(f,,903)  EPD 
VPITE(Gf90?)VPP 
PELr (ARP -EPP )/VRP**0.5 
WRITE  ( K,  9 I 9 )PEL. 


---ENP  OF  PAGE  43  — 


1U6 


1 


RUTPEP 


CONT 


I 0J  50|  40  04/21/78 


FILE  PAGE  NO, 


I I 


451  PC 

452  PC 
4530C 
454PC 
4550C 
4560 
4570 
45g0 
4590 
4600 
46 1 0 
4620 
4630 
464P 
465P 
466P 
467P 
4680 
469P 
47PP 
471  P 
472P 
4730 
4740 
475P 
4760 
4770 
47gp 
4790 
4B0P 
48IP 
48?P 
4830 
484  0 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
492P 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
503  0 
5040 
5050 


‘ --PAGE  *4 

ASPHAl TIC  CONCRETE  WITH  SUBGRAPF,  IP  VARIABLES. 


44  CONTINUE 

MP=MCI) 

MP=MP/I000. 

VP=V<1  ) 

VP=VP/|  000000. 

MTP =M(?) 

VTP=V(?) 

M?=M<3) 

VR=V<3) 

*TTI  =M(4) 

VT 1 =V(4 ) 

MT2=MC5> 

VT2rV(5) 

MCI =M(6) 

VCI =V(6) 

MT3=M(7) 

VT3=V(7) 

MC?=M(8) 

VC?  =V(  8) 

MC3=M(9) 

VC3:V(9) 

K =0.031 171 
A =1  .5255 
B =0.0897 
C =0.345 
P=-P.88^7 
E=-0 . 761 6 
F=- 1.1 674 
G=-0.5505 
H=-0.3089 

P 1 =F*  r»P** A*  MTP**B*  *R**C 
P?  = CAl  OG  I8<  1.25*  IT  I +MT?  >>**[) 

P3=MCI**  & (At  OG  10(MT3))**F*MC?**G*>C3**H 
P :*>  1*02*03 

P2=(A*CA-I)*VP)/MP**2 
TP2= (B* (B - 1)*  VTP ) /MTP**2 
P?  = <C* (C- 1 . )*  VR ) /MR**? 

T1?A=1.?5**?*0.434  3*P*VTI 

T 12B  =(  1 .25*  MT  1+  MT2  )**?*Al.OG10(  1 .25*  Ml  1 + IT?  ) 
TI2C  = ((0.a3a3*(D-1  ,))/AL0GI8(l  .25*  MT1  + MT2)) 
T1?=T12A/T 12B* (TI2C-I .) 

T??A=0.43A3*D*VT? 

T??=T??A/T 12B* (TI2C-1 .) 

C 1 2 = ( F>  (E-l  .)*VCI ) /MG  I **? 

T32A  =0  . A 34  3*  F*  VT  3 
T32B  =MT3**?*Al  OG  I0(MT3) 

T32C  =0.4  34  3*  (F-l  .J/At.OG  1 0(  MT3 ) 

T3?=T3?A /T32B* (T32C-1 ,) 

C2?  = <G*(G- 1 . )*  VC2 ) /MC?**2 
C3?= (H* (H- 1 . )* VC3 ) /MG3**2 


lUT 


I0»50J40 


04/21/78 


FILE  PARE  NO.  12 


RUTDEP  CONT 


9060  ERD:P+( 0.5*P*  (P2+TP2+P2+T1 2+  T22+C 1 2+  T32+C22+C32)  ) 

90  70  Plr(A**2*VP)/|*lP**2 

5000  TP  1 = <B**2*  VTP ) /MTP**2 

5090  P|  ::<C**2*W)/MP**2 

5100  Tl 1 A "I .25**2* <0. 4345* n>**2*VTl 

51  10  TUB  ;(1 .25*  Mri+'fr2)**2*(ALOG10(1.25*f|Tl  + MT2))**2 

5120  T!  1=T!  IA/TI  IP 

5130  T2  1 A ?( 0.4 34 3*0)**?*  VT2 

5140  T2  l=T2  1A0T1  IB 

5150  CM  =E**2*  VCI  /KCI**2 

5160  T31A  1(0.4 .54  J*  F)**2*  VT3 

5170  T3IP rMT3**2*Al  OG 10<MT3)**2 

5180  T3I:T31A/T3IP 

5190  C2I=(G**2*VC2)/MC2**2 

5200  C3 1 : (M**2*  VC3 )/MC3**2 

5210  VRD=P**2*  (Pl  + TP  l+PI  + TI  1 + T21+CI  1+T3  I+C2  1 -*-C  3 1 )-(0.25*P**2)* 

5220#,  (P2**2+TP2**2+P2**2+TI2**2+T?2**2+C  12**2+  T32**2+C22**2 

5230#,  -*-C32**2 ) 

5240  WPITE(6,904)P 

5250  WP I TF.  ( 6 , 9 03  ) *PD 

5260  WPITE(6,902)  VPD 

5270  PELs(APP-EPD)/VRD**0.5 

5280  WRITE (6, 9 19JPEL 

52  90C 

5300  RO  TO  1040 

53I0C 

5320  END 

5330C 

53  4 0C  --END  OF  PAGE  44--- 


RUTDFP 


CONT 


10»50j40  04/21/78 


FILE  PAGE  NO 


15 


53  5 00“ 
536PC 

53  7 PC 
538PC 
539P 
54PPC 

54  I PC 
54?  PC 
54  3 PC 
54  4 0 
54  5P 
54  6 0 
54  7PC 
548PC 

54  9 PC 
55PP 

55  IP 
55?  P 
553 P 15 
554p 
5550  10 
5560  800 
5570 
55gP  H 
5590 
56PP 
56IP 

56? 0 12 

563  PC 

564  PC 
565PC 
566P 
5670 
5680 

5690  180 
57PP 
5710 
57?  P 
5730 

5740  190 

575P 

576P 

5770 

5780 

5790C 

58PPC 

5SIPC 

5820 

5830 


— PAGE  SP- 
READ DATA  FROM  FILE.  COMPUTES  «EAN  A VAPIANCE. 
SUBROUTINE  STAT<  M,V,I ) 

READ  DATA  FROM  FILE.  COMPUTES  MEAN  A VARIANCE. 


DIMENSION  VALUE (300)  M( 1P),V(1P; 

INTEGER  TO  .GO 
PEAL  ?€AN,M 

INPUT  SCHEME  FOP  DATA  FILE. 

PEAD(IP,80O)LINENO,N 
GO:  1 
TO  =8 

CONTINUE 

IF  ( N-TO ) 1 P, 1 P, 1 1 

REAP (IP, 8PP) LI NENO, (VALUE (J),J=GO,N) 
FOPMAT(V) 

GO  TO  12 

PEAD  ( 10,800)  LI  NENO  , ( VALUE  (.1  > :GO  ,TO  ) 
GO=TO+l 
TO:TO+8 
GO  TO  15 
CONTINUE 

CALCULATE  fEAN  A VAPIANCE. 


SUMMrP. 

DO  IgP  J:1,N 

SUMM:VAl.UE(.l)+SUMM 
CONTI NUF 
MFAN:SUMM/N 
SUMVrP. 

DO  190  .1:1,  N 

SUMV=(VAl.UE(.l)-MEAN)**2+SUMV 
CONTI NUF 
IF ( N.GT .30) N:N- 1 
VAP  ^UMV/N 
M(I ):MFAN 
V(I)=VAP 

—END  OF  PAGE  50--- 

PFTUPN 

END 


lU9 


In  accordance  with  letter  from  DAEN-RDC , DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Barber,  Victor  C 

The  deterioration  and  reliability  of  pavements  / by  Victor  C. 
Barber,  Eugene  C.  Odom,  Robert  W.  Patrick.  Vicksburg,  Miss.  : 

U.  S.  Waterways  Experiment  Station  ; Springfield,  Va.  : 
available  from  National  Technical  Information  Service,  1978. 

xv,  149  p.  : ill.  ; 27  cm.  (Technical  report  - U.  S.  Army 
Engineer  Waterways  Experiment  Station  ; S-78-8) 

Prepared  for  Office,  Chief  of  Engineers,  U.  S.  Army,  Washing- 
ton, D.  C.,  and  U.  S.  Department  of  Agriculture  Forest  Service, 
Washington,  D.  C.,  under  Project  4A762719AT40,  Task  A2,  Work  Unit 
011,Q6,  and  USDA  Forest  Service-OCE  Interagency  Agreement. 
References:  p.  98-109. 

1.  Pavement  design.  2.  Pavement  deterioration.  3.  Pavements. 

4.  Reliability.  I.  Odom,  Eugene  C.,  joint  author.  II.  Patrick, 
Robert  W. , joint  author.  III.  United  States.  Army.  Corps  of 
Engineers.  IV.  United  States.  Forest  Service.  V.  Series:  United 
States.  Waterways  Experiment  Station,  Vicksburg,  Miss.  Technical 
report  ; S-78-8. 

TA7.W34  no. S-78-8 


